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1 Introduction 



The publish-or-perish period of science could soon be seen as a golden age: 
our brave new world now celebrates the triumph of Enron's financial and ac- 
counting creativity. Sadly science is now also, increasingly, considered from an 
accountant's viewpoint. In this respect, if one takes this "modern" point of view 
Professor F. Y. Wu's contribution^ is clearly a very good return on investment: 
he has given more than 270 talks in meetings or conferences, published over 200 
papers and monographs in rcfcrccd journals, and had many students. He has 
also published in, or is the editoij^ of, many books ^ [nj |l2|, |l3|, [l57, 



171] [L7§ ft79t [196 1 

Professor Wu was trained in theoretical condensed matter physics [|[ |[ [l9[ 



po| , [27], [35|, 108 1, but he is now seen as a mathematical physicist who is a leading 
expert in mathematical modeling of phase transition phenomena occurring in 
complex systems. Wu's research includes both theoretical studies and practical 
applications^. Among his recent researches he has studied connections of statis- 
tical mechanical models with deep problems in pure mathematics. This includes 
the generation of knot and link invariants from soluble models of statistical me- 
chanics and the study of the long-standing unsolved mathematical problem of 
multidimensional partitions of integers in number theory using a Potts model 
approach. 

Professor Wu's contributions to lattice statistical mechanics have been mostly 
in the area of exactly solvable lattice models. While integrable models have 
continued to occupy a prominent place in his work (such as the exact solution 
of two- and three-dimensional spin models and interacting dimer systems), his 
work has ranged over a wide variety of problems including exact lattice statistics 
in two and three dimensions, graph theory and combinatorics to mention just 
a few. His work in many-body theory |3|, |[ [?], ||, 15, p2[ |28|, 36 [S6|, especially 
those on liquid helium ||, ||, ^, ^6) , has also been influential for many years. 

F. Y. Wu joined the faculty of Northeastern University to work with Elliott 
Lieb in 1967, and in 1968 they published a joint papeiQ on the ground state 

1 Professor F. Y. Wu is presently the Matthews University Distinguished Professor at 
Northeastern University. He is a fellow of the American Physical Society and a permanent 
member of the Chinese Physical Society (Taipei). His research has been supported by the 
National Science Foundation since 1968, a rare accomplishment by itself in an environment of 
declining research support in the U.S., and he currently serves as the editor of three profes- 
sional journals: the Physica A, International Journal of Modern Physics B and the Modern 
Physics Letters B. 



2 For instance, Ref. [18C] contains the proceedings of the conference on "Exactly Soluble 
Models in Statistical Mechanics: Historical Perspectives and Current Status", held at North- 
eastern University in March 1996 - the first ever international conference to deal exclusively 
with this topic. The proceedings reflect the broad range of interest in exactly soluble models 
as well as the diverse fields in physics and mathematics that they connect. 

3 He has considered, for instance, the modeling of physical adsor ption an d applied it to 
describe processes used in chemical and environmental engineering [l48| , |l7E| ]. He has even 
published one experimental paper on slow neutron detectors |5|. 

4 This paper has become prominent in the theory of high-T c superconductors. P. W. 
Anderson even attributed to this paper as "predicting" the existence of quarks in his Physics 
Today (October, 1997) article on the centennial of the discovery of electrons. 
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of the Hubbard model [[LlJ which has since become a classic. The Baxter- Wu 
model fl9) is also, clearly, an important milestone in the history of integrable 
lattice models. 

F. Y. Wu has published several very important reviews of lattice statistical 
mechanics. First, Lieb and Wu wrote a monograph in 1970 on vertex models 
which has become the fundamental reference in the field for decades |3lJ . Wu's 
1982 review on the Potts model is another classic |8£j. At more than one 
hundred citations per year ever since it was published, it is one of the most 
cited papers in physics^]. In 1992 F.Y. Wu published yet another extremely 
well-received review on knot theory and its connection with lattice statistical 
mechanics [154|. In addition, in 1981, F. Y. Wu and Z. R. Yang published a 
series of expository papers on critical phenomena written in Chinese [ p4| - |88| . 
This review is well-known to Chinese researchers. 



1.1 The choice of presentation: A challenge in enumera- 
tive combinatorics 

An intriguing aspect of lattice statistics is that seemingly totally different prob- 
lems are sometimes related to each other, and that the solution of one problem 
can often lead to solving other outstanding unsolved problems. At first sight, 
most of the work of F. Y. Wu could be said to correspond to exact results 
in lattice statistical mechanics, but because of the relations between seemingly 
totally different problems it can cquivalently be seen, and sometimes be explic- 
itly presented, as exact results in various domains of mathematical physics or 
mathematics: sometimes exact results in graph theory, sometimes in enumera- 
tive combinatorics, sometimes in knot theory, sometimes in number theory, etc. 
Wu's "intellectual walk" goes from vertex models to circle theorems or duality 
relations, from dimers to Ising models and back, from percolations or animal 
problems to Potts models, from Potts models to the Whitney- Tutte Polynomi- 
als, to polychromatic polynomials or to knot theory, from results, or conjectures, 
on critical manifolds^ to Yang-Baxter integrability, perhaps on the way revis- 
iting duality or Lee- Yang zeros, etc., etc. The simple listing of Professor Wu's 
results and contributions, and the inter-relations between these results and the 
associated concepts and tools, is by itself a challenge in enumerative combina- 
torics. 

Actually it is impossible to describe Wu's contributions linearly, in a se- 
quence of sections in a review paper like this, or even with a website-like "tree 
organization" of paragraphs. F. Y. Wu's contributions really correspond to a 
quite large "graph" of concepts, results, tools and models, with many "intel- 
lectual loops" . The only possible "linear" and exhaustive description of Wu's 
contributions is his list of publications. 

5 There was once a study published in 1984 (E. Garfield, Current Comments 48, 3 (1984)) 
on citations in physics for the year of 1982. It reports that in 1982, the year this Potts review 
was published, it was the fifth most-cited paper among papers published in all physics. 

6 The critical manifolds deduced or conjectured by F. Y. Wu are mostly algebraic varieties 
and not simple differentiable or analytical manifolds. 



5 



We have therefore chosen to give his exhaustive list of publications at the end 
of this paper. No other references are given. 

We have chosen to keep the notation F.Y. Wu used in his publications^ and 
not to normalize them, so that the reader who wants to see more and goes back 
to the cited publications will immediately recover the equations and notations. 

Obviously, wc will not try to provide an exhaustive description of Wu's con- 
tributions but, rather, to provide some considered well-suited specific "morceaux 
choisis^" , comments on some of his results, some hints of the kind of concepts 
he likes to work with, and trying to explain why his results are important, fruit- 
ful and stimulating for anyone who works in lattice statistical mechanics or in 
mathematical physics. 

2 Even before vertex models: the exact solution 
of the Hubbard model 

Elliott H. Lieb and F. Y. Wu published in 1968 a joint paper on the ground 
state of the Hubbard model jll) which has since become a classic, and served as 
a cornerstone in the theory of high-T c superconductors. An important question 
there corresponds to the spin-charge decoupling which is exact and explicit in 
one-dimensional models: is the spin-charge decoupling a characteristic of one 
dimension? Is it possible that some "trace" of spin-charge decoupling remains 
for quantum two-dimensional models which are supposedly related to high-T" c 
superconductors? 

Let us describe briefly the classic Lieb-Wu solution of the Hubbard model. 
One assumes that the electrons can hop between the Wannier states of neigh- 
boring lattice sites and that each site is capable of accommodating two electrons 
of opposite spins with an interaction energy U > 0. The corresponding Hamil- 
tonian reads: 

H = T Yl H c ^ + u Yl c U Ci t c U c u 

<ij> ct i 

where c\ a and Cj „ are the creation and annihilation operators for an electron 
of spin a in the Wannier state at the i — th lattice site and the first sum is 
taken over nearest neighbor sites. Denoting f(xi, x%, • • • , xm\ £a/+i, • • • ; %n) 
the amplitude of the wavefunction for which the down spins are located at 

7 The price paid is, for instance, that the spin edge Boltzmann weights will sometimes be 
denoted e Kl , e K2 , e K z , e K 4 , or a, b, c, d, or xi , X2,%3, £4, and the vertex Boltzmann weights 
u>i, a)?, ■ ■ ■ or o, b, c, d, a' , 6', c', d' . This corresponds to the spectrum of notations used in the 
lattice statistical mechanics literature. These different notations were often introduced when 
one faces large polynomial expressions and that the e Ki or e~@' Ji notations for Boltzmann 
weights would be painful. 

8 I apologize, in advance, for the fact that these "morceaux choisis" are obviously biased by 
my personal taste for effective birational algebraic geometry in lattice statistical mechanics. 
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sites xi, X2, ■ ■ ■ ,Xm and the up spins are located at sites xm+i, ■ ■ ■ ,£at- The 
eigenvalue equation Hip — E ip leads to: 

N 

-^2^2 f{x!, X 2 , ■ ■ ■ ,Xi + S, ■ ■ ■ ,X N ) (1) 
i=l s =±l 

+ J7 } S(xj -Xj)f(xi, x 2 , ■■■ ,x N ) = Ef(xi, x 2 , ■■■ ,x N ) 

where f{x±, x 2 , ■ ■ • ,xn) is antisymmetric in the first M and the last N — M 
variables (separately). Let /x+ (resp. denote the chemical potential of 

adding (resp. removing) one electron. In the half-filled band one has /i+ = 
U — /Lt_, and the calculation of //_ can be done in closed form with the result: 

n =2 _ 4 MUj) ' dLJ (2) 

' Jo u ■ (1 + exp (w C//2)) v y 

where Ji is the Bessel function. It can be established from (j^) and /i + = 
J7 — /i_ that //_|_ > /i_ for [/ > 0. In other words, the ground state for a 
half-filled band is insulating for any nonzero U, and conducting for U = 0. 
Equivalently, there is no Mott transition for nonzero U, i.e., the ground state 
is analytic in U on the real axis except at the origin. 



3 Vertex models 

The distinction between vertex models and spin models is traditional in lattice 
statistical mechanics, but there are "bridges" between these two sets of lattice 
models Q. Roughly speaking one can say that F. Y. Wu first obtained re- 
sults on vertex models Jl3|, [l4| (five- vertex models ||, |l(J, free-fermion vertex 
models |50[| , dimer models seen as vertex models, ...) and then obtained re- 
sults on spin models (Ising model with second-neighbor Interactions |Q, the 



Baxter-Wu model 45, |49|, Potts model, ...), introducing more and more graph 
theoretical approaches up to looping the loop with knot theory which is, in fact, 
closely related to vertex models and Potts models! As far as vertex models are 
concerned, we will first sketch the approach given in his monograph with Lieb 



(section (3.1)), in a second step we will sketch his free-fermion results (section 
( [3.2.1 )) closely followed by his dimer results (section (3.3)), and, then, we will 



discuss some miscellaneous results he obtained on five-, six- and eight- vertex 



models (section (3.4)) 



3.1 Two-dimensional Ferroelectric Models 

Elliott Lieb and F. Y. Wu wrote a monograph on vertex models in 1970, entitled 
"Two-dimensional Ferroelectric Models" , which has become a fundamental ref- 
erence in the field for decades [|l| . This monograph gives the best introduction 
to the sixteen-vertex model, which is a fundamental model in lattice statistical 
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mechanics. But unfortunately it is not known well enough even to many spe- 
cialists of lattice models that it contains the most general eight-vertex model, 
most of the (Yang-Baxter) integrable vertex models (the symmetric eight- vertex 
model, various free-fermion models, the asymmetric free-fermion model, the 
asymmetric and symmetric six-vertex, the five-vertex models, three-coloring 
of square maps, and others) and also fundamental non-integrable models such 
as, for instance, the Ising model in a magnetic field. In particular the mono- 
graph mentions explicitly the weak-graph duality (see section (|^) below) on the 
sixteen- vertex model (see page 457 of Q): 



The 154 pages of this monograph are still, by today's standard, an extremely 
valuable document for any specialist of lattice models. Beyond the taxonomy 
of ferro and ferrielectric models (ice model, KDP || |l8), modified KDP ifflf , 
F model modified F model (3^, [75|, p(J, F model with a staggered field, 
...), this monograph remains extremely modern and valuable from a technical 
viewpoint. 

Among the exactly soluble models (the bread-and-butter of F. Y. Wu) was 
one that, for a long time, was a "sleeper", namely, Bethe's 1931 solution of 
the ground state energy and elementary excitations of the one-dimensional 
quantum-mechanical spin-i Heisenberg model of antifcrromagnetism. We will 
see below a large set of results from the Lieb-Wu monograph on vertex models, 
in particular the six-vertex model. 

The monograph gives an extremely lucid exposition of the Bethe ansatz for 
the six-vertex model. The Bethc ansatz is analyzed and explained in the most 
general framework (with horizontal and vertical fields) and it is a must-read 
to anyone who wants to work seriously on the coordinate Bethe ansatz. It is 
certainly much more interesting and deeper than so many subsequent papers 
that have revisited, at nauseum, the Bethe ansatz of the symmetric six- vertex 
model, re-styling this simple Bethe ansatz with a conformal resp. quantum 
group, resp. knot theory, resp. ... framework. The analysis of the conditions 
for the transfer matrix T of the most general sixteen-vertex model to have a 
non-trivial "linear operator" (ID quantum Hamiltonian) that commutes^ with 
T (pages 367 to 373) are probably one of the first pages any student who wants 
to study integrable lattice models should read. 

The monograph makes crystal clear of the fact that the Bethe ansatz is 
related to the conservation of certain charge. This is seen from the fact that 
most of the analyses (from page 374 to page 444) relies on the use (page 363 

9 Which is the most obvious manifestation of the Yang-Baxter integrability. 
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equation (81)) of the variable y = 1 — 2 n/N, which in the spin language is the 
average z-component of the spin per vertical bond, namely, y = < S z > /N for 
a square lattice of size N X M , where n denotes the number of down arrows 
and iV the number of vertical bonds in a row. 

We use the same notation as in Lieb-Wu. In particular let us introduce the 
horizontal and vertical fields H and V respectively. The partition function per 
site in the thermodynamic limit is: 

Urn — • ln(A) = max \z(y) + V ■ y] (4) 

N —* oo iv — l<y<+l 

where A denotes the largest eigenvalue of the transfer matrix. The monograph 
details a large set of situations. Let us consider here the regime 

A = (wiw 2 +W3W4 — cu 5 uj 6 ) /2-y/wi 0J2 ^3 ^4 < — 1, (5) 

and introduce the variable: 

e e ° = \ + 77 gA , < < A where: 



= e \ Kl \ = f^i^y/ 2 or . /^4_y/ 2 

When A < — 1, z(y) reads: 

1 f +b 

z(y) = -K 2 + max(0, -Kx) ^ / • C(a) • da 

- n f \ 1 /cosh(2 X — 9 ) — cos(a)\ 

where: C{a) = Znl 

V cosh(t'o) — cos (a) / 

and the (normalized) density^] R(a) satisfies the Bethe-ansatz integral equation 
with the kernel K(a); 

R(a) = — u . S x mh(A:> - - / + K(a - /3) ■ R{fi) dp (6) 



cosh(A) — cos(o;) 
with: 2 7T • K(a - 0) 



-b 

sinh(2 A) 



cosh(2 A) — cos(a — 0) 



The integral equation (||) is nothing but the well-known Yang- Yang Bethe ansatz 
integral equation on the density p(q) : 

Q dd(p,q) , v j • 1 1 ^ TT.(l-tf) 



1 = 2 7T • p(p) - / ^M±p( q ) . dq with: Q 

7-q dp 



The range b of the new variable a in the integral relation (^) can be deduced 
from the definition of the density R(a) : 

r+b r+Q 

7T • (1 — y) = I R(a) da = / p(q) ■ dq . 

J-b J-Q 

10 Ono has R(a) ■ da = 2 n p(p) ■ dp. 



When y = the integral attains its maximum range and one can solve (||) by 
a Fourier series of a Fourier transform. One thus gets R(a) as a simple dn 
elliptic function. Not surprisingly one can also calculate all the derivatives of 
z(y) at y = . One can thus expand z(y) namely, write z(y) = z(0) — z'(0)- 
y + z (0) • j/ 3 /6 + • • • . In first order in y one obtains z'(0) = — E(A — 9q) 
where the function S is related to the Jacobian elliptic function nd: 



_ / COsh((A + 6)/2) \ _ £ _ y, (_l)". e -2nA , smh(n 

" W HcoshffA - 6)12)) 2 ^ n-cosh(nA) 



-2 n\ 

sh((X - 4>)/2)J 2 ^ n ■ cosh(n A) 

The function H(^) also satisfies the nice involutive functional relations^] 

E(6) = -E(-6), 3(A + 0) = 5(A - 6), S(4 A + 6) = E(6). 

Let us consider the thermodynamic properties of the model when H = 
and V 7^ 0. From (Q) one sees that the thermodynamic properties depend on 
the optimal choice of y given by: 

z'(y) = - V 

When lowering the temperature the slope of z{y) corresponding to the tran- 
sition sticks at y ~ 0, and one thus has (see page 425 of J3l|]) an antiferroelectric 
transition occurring at T C (V) given by: 

V = 5(A - 9 Q ). (7) 

This gives a beautiful example of a transcendental critical manifold which 
reduces, in some domain of the parameters (low temperatures), to a transcen- 
dental equation (Q) and not to an algebraic one as one is used to see in exactly 
solvable models. One thus has a transcendental critical manifold for a vertex 
model for which one can actually write down the exact Bethe ansatz (see equa- 
tion (||)). Writing a closed simple formula for the solution is not possible, but 
one can certainly find numerical solutions on a computer. Should we say that 
the model is exactly solvable but not "computable" ? We will revisit these ques- 
tions of algebraicity of the critical manifold versus integrability in other sections 
of this paper with other criti ca l ma nifolds conjectures, or results, of F. Y. Wu 
(see for instance sections (3.3), ( |6.l[ ) below). For those who have a "naive" point 



of view on the character of critical manifolds^], the example (0) shows that a 
model having a Bethe ansatz can have a transcendental critical manifold. 

The Lieb-Wu review provides wonderful pieces of analytical work (analysis 
in one complex variable, see for instance pages 410-411 and the analysis of the 

11 In agreement with the inversion relations on the model. 

12 With, for instance, a prejudice of algebraicity of the critical manifolds of "solvable" models: 
all examples known in the literature are polynomial expressions in well-suited variables e Ki . 
These include, for instance, the critical varieties of the anisotropic Ising, or Potts, models on 
square, triangular lattices, or the critical varieties of the Baxter model. For non-integrable 
models the common wisdom is, probably, that critical manifolds are always analytic, or may be 
differentiable, and the algebraicity of the critical manif olds is ruled out by the non-integrability. 



This is also a naive point of view: see (|36|) in section (5.1) 
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analytic structure of the F model or the temperature Riemann structure for the 
free energy of the F model). One finds a festival of one complex variable analyt- 
ical tools (Maclaurin formula, tools for the evaluations of asymptotic behaviors, 
path integration, etc.). 

Many more results can be found in the monograph (the three-color problem, 
the hard square model, the F model on the triangular lattice, three coloring 
of the edges of the hexagonal lattice ...). Let us mention, in particular, the 
six-vertex model with site- dependent weights (which can be considered as the 
first example of a Z-invariant model). Let us introduce u>j(I,J) where j = 
1, • • • , 6, are the six possible Boltzmann factors of the vertex in row / and 
column J, and let us require that the algebraic invariant A be independent of 
I, J: 

a = m(/,JWJ,J) + - u 5 (i,J)uj 6 (i,J) 

2-{u 1 {I,J)u 2 {I,J)uz(I,J)u H {I,J))V* ■ U 

Up to a multiplicative factor Kj t j , a (rational) parametrization of these invari- 
ance conditions (g) is: 

Ui(I,J) = (1 - t -pi,j) • a/,j/3/,j, 

lj 2 (I, J) = (1 - t ■ Pi.j) (9) 

di, j pi,j 

u 3 (I,J) = {pi,j - t) ■ w 4 (J, J) = (p/,j - t) ■ 

PI,J OLI,J 

w B (/, J) = (t - t) ■ pi.j ■ 7 /,J, w 6 (I, J) = (1 — t 2 ) • — 

* H,J 

The Baxter's Z-invariance condition for integrability requires that the pi,fs are 
actually products of a (spectral) parameter depending on the row and another 
parameter depending on the column: pj t j = pj ■ oj. We will see in section (0) 
when sketching the correspondence between the standard scalar Potts model 
and a staggered asymmetric six- vertex model, that these product conditions, 
Px j = pj <tj, actually correspond in the case of the checkerboard Potts model 
to criticality, or to a vanishing conditions of a staggering field H stlig in a Lee- 
Yang zeros viewpoint: \z\ = e ffstag = 1. 

Provided that the pi t j = pj ■ aj integrability conditions are satisfied, the 
partition function, with parametrization (^) , can be expressed as a multiplicative 
closed formula: 



i=i i=i 



where: F(z) = JJ j—^+r^ (10) 

m. — 1 



n 



3.2 Vertex models: free fermions 

Another classic work of F. Y. Wu is his 1970 paper with C. Fan in which they 
coined the term the free-fermion model ]l6[ . This work was later extended 
to its checkerboard version during one of Wu's visit to Taiwan |5(| |5^]. In the 
following we shall arrange the homogeneous vertex weights in a matrix R, whose 
size and form vary according to the number of edge states and the coordination 
number of the lattice. Typical examples we will consider are the 2D square and 
triangular lattices shown below: 



v 

W V 




3 3 k 

2D square 2D triangular 



3.2.1 Free-fermion asymmetric eight-vertex model 

C. Fan and F. Y. Wu obtained many free-fermion results [[[2], The free 

energy of the most general free-fermion model on a square lattice evaluated by 
Fan and Wu reads: 



1 



16 tH 



(i 



2, 

de d</> In ( 2 a + 2 b cos6> + 
2 c cos 6 + 2 d cos(9 — 6) + 2 e cos( 



where: 



I 2 i 2 , 2 , 2\ ; 

a = - ■ [uj-y + lu 2 + oj 3 + U! 4 ), b = uji uj 3 — w 2 <^i, 

C = Ll>i LU4 — L02 U>3, d — LU3 UJ4 — UJj Wg , € = LU3 LU4 — U>Q , 



provided that the free fermion condition: 

Wl ^2 + UJ3UJ4 = W 5 W 6 + W7W8 (11) 

is satisfied. 

Let us revisit some of their results from an inversion relation viewpoint. 
Renaming the vertex weights as a = u>\, a' = u>2, b = lo 3 , b' = lo 4 , c = u> 5 , 
d = luq, d = lut, d' — uig, the matrix R of the eight-vertex model is then: 

fa cT 
\d d 
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A matrix of the form ( |12| ) can be brought, by a similarity transformation, to a 
block-diagonal form: 

*-(* I), wit, *,= (« *) and K 2 = (^ =!) 

If one introduces <5i = aa! — dd' and ^ = bb' — cc', the determinants of the two 
blocks, then the (homogeneous) matrix inverse I (namely R — > det(i?) • R^ 1 ) 
reads: 

{a, a! , d, d') — > (a' • £2, a • 82, — d- 82, — d! ■ 82) (13) 
(6, b', c, c') -» (6'-5i, &-<5i, -c-5i, -c'-<5i) 

It is straightforward to see that the free-fermion condition (^lj) is <5i = —82 
which has the effect of linearizing the inversion ( |l3| ) into an involution given by: 

a^a', b^-b', (d,d')~>(-d,d'), b^-b', {c, c') ^ (c, c') 

The group generated by the two inversion relations of the model is then realized 
by permutations of the entries mixing with sign changes, and its orbits are 
thus finite. The finiteness condition of the group is a common feature of all 
free-fermion models. 



3.2.2 Free-fermion for the 32-vertex model on a triangular lattice 



We next consider the free-fermion conditions of J. E. Sacco and F. Y. Wu |33| 
for the 32-vertex model on a triangular lattice. Using the same notation as 
in we have: 



R 



fo 








}23 





fl3 


fl2 








fee 


f-2 6 





fl6 








fj.5 





fs5 


f25 





fl5 








fl,C> 


f56 








fu 





fz4 


f34 








fs4 


f24 





fu 








f56 










fit 5 





fj25 


fs5 













fl6 





f26 


fs6 








fl2 


fl3 





}23 








fo 



(14) 



By permuting rows and columns, this matrix can be brought into the block 
diagonal form: 



R = 



i?i 
R 2 



with: 



Ri = 



fo 


ft 3 


fl2 


f23 


f46 


$25 


f§5 


flB 


f45 


f&6 


f§6 


fl6 


fs6 


fz4 


f34 


fl4 



Ri = 



fu, 


fs4 


f&4 


fee 


fl6 


ha 


f&6 


fj5 


f_15 


fs5 


f&5 


f46 


f23 


}l2 


}l3 


fo 



(15) 



(16) 



The inverse / written polynomially (homogeneous matrix inverse) is now a 
transformation of degree 7. If one introduces the two determinants Ai = 
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det(i?x) and A 2 = det(i? 2 ), then each term in the expression of I(R) is a 
product of a degree three minor, taken within a block, times the determinant 
of the other block. This inverse I clearly singles out one of the three directions 
of the triangular lattice. These three involutions do not commute and generate 
a quite large infinite discrete group r tr i a ng (see also section ( |6.l|) below). 
The free-fermion conditions of Sacco and Wu |53| read: 

fofijkl — fijfkl — fikfjl + filfjk, V i, j, k, I = 1, . . . , 6 

/o/o = /12/12 — /13/13 + /14/14 — /15/15 + /16./16 (17) 

which we denote by V. What is remarkable is that, not only is the rational va- 
riety V globally invariant under r tr iang, but again the realization of this (gener- 
ically very large infinite discrete group) r tr i a ng on this variety becomes finite. 
This comes about from the degeneration of / into a mixture of sign changes 
and permutations of the entries, as in the case in the preceding subsection. 

Remark: The ordinary matrix product of three matrices ( |l4] ) solutions of ( [i"7| ) 
is another solution ! In other words, if R a , Rp, i? 7 € V, then R a ■ Rp ■ i? 7 6 V, 
while R a ■ Rp V. This was also the case for involutions of ( |TT| ) in the case 
of the square lattice, but the mechanism is more subtle here as conditions ( p~7| ) 
imply Ai = A 2 . 



3.3 Dimers and spanning trees 

Before Fan and Wu's free-fermion vertex models the Onsager solution of the 
two-dimensional Ising model is clearly the first free-fermion model ever solved. 
There were also several approaches to the two-dimensional Ising model that 
did not use the transfer matrix formalism; the most interesting one is perhaps 
the mapping of the problem onto a dimer-covering problem on a slightly more 
complicated lattice. The dimer problem was first solved by Temperley-Fisher 
and Kasteleyn. Kasteleyn found out how to treat the most general planar graph. 

The dimer problem has a life of its own and has since found many followup 
works, not only in statistical mechanics, but also in combinatorial theory. In 



this rega rd Wu has provided a large number of new results [B71 173, 184, 194 



207 , 205 j, including applications to condensed matter physics as well as in pure 
combinatorial anal ysis . In a ddition, Wu has obtained new results on the span- 
ning tree problem |l98| , [200| , a problem intimately related to the dimer problem 
through a bijection due to Temperley. In the following we shall describe some 
of the contributions in this area. 



3.3.1 Revisiting dimers: the honeycomb lattice 

The dimer model on the honeycomb lattice was first solved by Kasteleyn, but he 
never published the solution except hinting at the existence of a transition. This 
deficiency was made up by Wu in a 1968 paper j[o) in which he presented details 
of the analysis for the honeycomb lattice, and applied the results to describe 
the physics of a modified KDP model. 
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3.3.2 Revisiting dimers: Interacting dimers in 2 and 3 dimensions 



Almost 30 years after the publication of the solution of the dimers on the hon- 
eycomb lattice [|l0] , Wu and co-workers made two important extensions of the 
earlier Kasteleyn solution. In the first, H. Y. Huang, F. Y. Wu, H. Kunz, and 
D. Kim | |173fl considered the case where the dimers have nearest-neighbor in- 
teraction. This model turns out to be identical to the most general five-vertex 
model, a degenerate case of the six- vertex model which requires a special Bethe 
ansatz analysis. The resulting phase diagram of this five-vertex model is very 
complicated and the analysis extremely lengthy. 



In the second work H. Y. Huang, V. Popkov and F. Y. Wu 177 , 184] intro- 
duced, and solved, a three-dimensional model consisting of layered honeycomb 
dimer lattices as described in the preceding subsection, but with a specific layer- 
layer interaction. Again, the phase diagram is very complicated. It is noted 
that this model is the only solvable three-dimensional lattice model with phys- 
ical Boltzmann weights (the Baxter solution of the 3D Zamolodchikov model 
has negative weights). However, the layered dimer model, while having strictly 
positive weights, describes dimer configurations in which dimers are confined in 
planes. As a consequence the critical behavior is essentially two-dimensional. 



3.3.3 Revisiting dimers: a continuous-line model 



F. Y. Wu and H. Y. Huang [ 158 have further used a dimer mapping to solve 
a continuous-line lattice model in three and hig her d imensions. They have also 
applied it to model a type-II superconductor [ 160 1 . In three dimensions the 
model is a special case of an 0(n) model on a finite L\ x L2 x L3 cubic lattice 
with periodic boundary conditions with the partition function: 



Z(n) 



E 

closed polygons 



n l ■ z\ 



where the summation is taken over all closed non-intersecting polygonal config- 
urations, I is the number of polygons, and b is the number of edges of each 
configuration. They considered the n = — 1 special case which they showed to 
be in one-to-one correspondence with a dimer problem whose partition function 
can be evaluated as a Pfaffian. The result for a finite lattices is: 



Z(-l) 



L\ L2 L3 

n n n 

m —1 712 —1 ™3 —1 



i=l 



In the thermodynamic limit, this leads to the per-site free energy: 



/ 



1 



2 TT 



2 TT 



2 TT 



dO* In 



(2tt)3 jo jo Jo 

The phase diagram is rich and quite non-trivial. 



3 

E 

i=l 



_ 2 tt m/Ni 
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However, it must be said that this exactly soluble three-dimensional 0(— 1) 
model describes line configurations running only in a preferred directions and, 
secondly,^ the Boltzmann weights can be negative. 



3.3.4 Revisiting dimers: nonorientable surfaces 



More recently W. T. Lu and F. Y. Wu initiat ed studies on dimers and Ising 
models on nonorientable surfaces [194, [203 , 20E ]. For dimers on an M x N net, 
embedded on non-orientable surfaces, they solved both the Mobius strip and the 
Klein bottle problems for all sizes M and N and obtained the dimer generating 



function Z 



M,N 



Z M N/2 

= Re(( 



(M+l)/2 AT 

n n( 2< (- i ) Af/a+ * n+isiii ( 



(4 n — 1) 7r 
2N 



2 X „ 



where Re denotes the real part, z v and Zh are the dimer weights in the vertical 
and horizontal directions respectively, and X m is given for the Mobius strip and 
the Klein bottle respectively by: 



z v ( m 7T 

X m = — - cost 

Zh \ M + 1 



z v /(2m- l)?r 
X m = — ■ cos( — 

z h \ M 



In paper [205] they also obtained an extension of the Stanley- Propp reciprocity 
theorem for dimersp*]. Inspired by this work, there is now much activity in this 
area. There is also very much current interest in finite-size corrections and 
conformal field theories on more complicated surfaces (higher genus pretzels, 
...). 



3.3.5 Dimers on a square lattice with a boundary defect 

In a very recent paper [ p07[ fittingly dedicated to the 70th birthday of Michael 
Fisher who first solved the dimer problem for the square lattice, W. J. Tzcng 
and F. Y. Wu obtained the dimer generating function for the square lattice with 
one corner (or some other boundary site) of the lattice missing. In this work 
they made use of a bijection between dimer and spanning tree configurations 
due to Temperley (and extended by Wu in his unpublished 1976 lecture notes 
as well as more recently by Kenyon, Propp and Wilson). They also carried out 
finite-size analyses which lead to a logarithmic correction term in the large-size 
expansion for the vacancy problem with free boundary conditions. They found a 
central charge c = —2 for the vacancy problem, to be compared with c = — 1 
when there is no vacancy. This central charge c = —2 is in contradiction with 
the prediction of the conformal field theory. 

13 But we are used to this after R. J. Baxter's solution of the 3D Zamolodchikov model. 
14 A subject matter of pertinent interest to mathematicians. 
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3.3.6 Spanning trees 



As mentioned above, the problem of spanning trees in graph theory is intimately 
related to the dimer problem, and it is not surprising that Wu found his way to 
spanning trees. In 1977 he published a paper J62[ on the counting of spanning 
trees on two-dimensional lattices using the equivalence with a Potts model. 
Very recently he refined the tools by using a result in algebraic graph theory 
which he and W. J. Tzeng rederived using elementary means. Tzeng and Wu 
enumerated spanni ng t rees for general d-dimensional lattices as well as non- 
orientable surfaces [ 198 1 . Applyin g th ese results to general graphs and regular 
lattices, R. Shrock and F. Y. Wu |20C|] published a lengthy paper in which they 
established new theorems on spanning trees as well as enumerating spanning 
trees for a large number of regular lattices in the thermodynamic limit. 



3.4 Miscellaneous results on vertex models 

In this section we describe an arbitrary choice of miscellaneous results obtained 
by F. Y. Wu on vertex models. 



3.4.1 Boundary conditions 

The six-vertex model is known to be a boundary condition dependent model. 
However H. J. Brascamp, H. Kunz and F. Y. Wu Ej| have established, for the 
first time, that, at sufficiently low temperatures or sufficiently high fields, the 
six- vertex model with either periodic or free boundary conditions are equivalent. 



3.4.2 The eight-vertex model in a field 



A simple result due to F. Y. Wu [ 105 1 is that a very general staggered eight- 
vertex model in the Ising language (as introduced by Kadanoff and Wegner and 
by F. Y. Wu |7£|) with the special Yang-Lee magnetic field iirkT/2, is equivalent 
to Baxter's symmetric eight-vertex model and hence is soluble. This result is 
remarkable since the general eight-vertex model without this field is not known 
to be soluble. 



3.4.3 The eight-vertex model on the honeycomb lattice 



F. Y . Wu has always been keen in providing results for the honeycomb lattice J48], 
130 1 . One interesting result is that he has established the exact equivalence 



of the eight-vertex model on the honeycomb lattice with an Ising model in a 
nonzero magnetic field M, 13C ] . The equivalence also leads to exact analysis 



of the Blume-Emery-Griffiths model for the honeycomb lattice (for details see 
section (4.3) below). 

In pursuit of applications of these results, P. Pant and J. H. Barry and 



Wu [181, 182 1 obtained exact results for a model of a ternary polymer mixture 
which is equivalent to an eight-vertex model. A model ternary polymer mixture 
was considered with bi- and tri-functional monomers and a solvent placed on the 
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sites of a honeycomb lattice. Using the equivalence with an eight-vertex model 
which further maps the problem into an Ising model in a nonzero magnetic 
field, an exact analysis of the model was carried out. The phase boundary of 
the three-phase equilibrium polymerization regime was determined exactly. 

Comment: These kinds of results are particularly interesting when one real- 
izes that concepts and structures corresponding to the Yang-Baxter integrability 
do not exist. This leads to the natural question: how to construct a Yang- Baxter 
relation for the honeycomb lattice? 



3.4.4 Exact critical line of a vertex model in 3 dimensions 

Wu |Jl6f has introduced a vertex model in three dimensions with real vertex 
weights, and determined its exact first-order phase transition line by mapping 
it to an Ising model in a field. It also exhibits a critical point. This is one of 
the very few lattice statistical models for which exact results can be deduced in 
higher-than-two dimensions. 



4 Spin models: Ising models and other models 

We will consider in this section F. Y. Wu's results on spin models, mostly Ising 
models. Due to its importance the Potts models will be treated separately in 
section ([?]). 

The distinction between vertex models and spin models (more generally 
Interaction Round a Face (IRF) models) is an important one in lattice statistical 
mechanics. However Wu showed in paper |7q] an equivalence between an Ising 



model with a vertex model. Similarly in paper \114\ Wu and K. Y. Lin studied 



the Ising model on the Union Jack lattice showing it to be a frcc-fcrmion model. 



Many of the free-fermion results on vertex models of sections ( 3.2.1 ) and (3.3) 
can also be re-styled as free-fermion Ising models. 

As far as the Union Jack lattice is concerned Wu has also obtained the spon- 
taneous magnetization of the three-spin Ising model [filf . It is, in fact, obtained 
in terms of the magnetic and ferroelectric orderings of the eight-vertex model, 
or, equivalently, the spontaneous magnetization and polarization of the eight- 
vertex model. It was found that the two sublattices possess different critical 
exponents. 

An important development in the history of lattice models is the analysis of 
the phase diagram of the Ashkin- Teller model on the square lattice by F. Y. Wu 
and K. Y. Lin ^7j. The Ashkin- Teller model is another example of spin models 
for which the traditional distinction of lattice statistical mechanics between spin 
and vertex models is irrelevant. The Ashkin- Teller model can be seen as two 
Ising models coupled together with four-spin interactions. Performing a dual 
transformation on one of the two Ising models and interpreting the result as a 
vertex model one finds that the Ashkin- Teller model is equivalent to a staggered 
eight-vertex model p9[, thus exhibiting two phase transitions. 
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Wu's analysis of spin models are not restricted to two dimensions. For 
instance, Barry and Wu have obtained exact results for a four-spin-interaction 



Ising model on the three-dimensional pyrochlore lattice 12S], and Wu has also 
obtained various results for spin models on Bethe lattice and Cayley trees jsi], 



Wu also performed real-space renormalization studies for Ising models [ 129 1 , 
but, not surprisingly, using some duality ideas, namely, the duality-decimation 
transformation of T. W. Burkhardt. Wu had previously applied the duality- 
decimation transformation in order to solve the two-dimensional Ising model 
with nearest-neighbor, next-nearcst-ncighbor and four-spin interactions in a 



pure imaginary field |l0§ (see section . In paper [|129| Burkhardt 's method 



which combines a bond-moving and duality-decimation transformation is mod- 
ified in order to preserve the free energy in the renormalization transformation. 

4.1 Generalized transmissivities for spin models 

We will see below that a large number of F. Y. Wu's work correspond to graph 
expansions (see section ( ]7.1j )). For spin models with edge interactions this re- 
quires the introduction of certain "transmissivity" variables. Thermal trans- 
missivities are introduced when considering high-temperature expansions of an 
edge-interaction spin model or performing renormalization analyses. They are 
also the natural variables to use in the decimation of spins in a simple multi- 
plicative way. 

Introduce the edge Boltzmann weight W(K±, K2, ■ ■ ■ K n ; a, b) where K\, 
K2, ■ ■ ■ K n denote a set of coupling constants describing the model, and a and 
b two nearest-neighbor spin states which can take on q values. Let us assume 
that the decimation procedure yields a Boltzmann weight of the same form: 

Y,W{K 1 ,K 2 ,---K n ; a, b) ■ W(K[, K' 2 , ■ ■ ■ K' n , 6, c) 

b 

= \-W{K'{, K'l a, c) (18) 

Alternatively, one can build a q x q Boltzmann matrix W with entries W^j = 
W(Ki, K2, ■ ■ ■ K n ; i, j). In terms of such matrices the relation Jl8| ) becomes 
W • W = W". The decimation procedures and also the high-temperature ex- 
pansions in such models are greatly simplified by introducing a "transmissivity" 
function t a such that the matrix relation W • W = W" becomes one or more 
multiplicative relations of the form: 

t a (W)-t a (W) = t a (W"), a = l,--.,r. 

The simplest example is the transmissivity variable for the q state standard 
scalar Potts model, for which one has t — (e K — l)/(e A + q — 1). This 
is the natural expansion variable for the high-temperature s eries of the model 



(see also the /i,j's in (|45| ) and (|46| ) introduced in section (7T) below). For 



the Ising model this reduces to the tanh(isr) variable. F. Y. Wu et al. [ 147 1 



underlined the fact that two quite different situations must be considered. If 
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the family of Boltzmann matrices W is a set of commuting matrices, then 
they can be diagonalized simultaneously and the transmissivity variables are 
nothing but all the possible ratios of eigenvalues of the Boltzmann matrices 
W. If, alternatively, the Boltzmann matrices W do not commute, then one 
must perform a simultaneous block-diagonalization of this family of Boltzmann 
matrices, and, therefore, some of the t Q 's will be block matrices from which 
one can extract functions 4> a satisfying 4> a (W") = <f> a (W) ■ <f> a (W'). One 
obvious choice for <p a is the (ratio of) determinants of these blocks. A number 
of non-trivial non-commuting transmissivities are given in [ 147 1 . 



4.2 Three-spin interactions: the Baxter- Wu model 

Another important work in the history of exact solutions of lattice statistics is 
the Baxter- Wu model, which is an Ising model on the triangular lattice with 
three-spin interactions. This model was solved exactly by R. J. Baxter and F. 
Y. Wu in 1973 

The three spin surrounding every triangular face interact with a three-body 
interaction of strength — J, so that the Hamiltonian reads: 

H = -J-J2<?i 03 °k (19) 

Baxter and Wu found that the per-site partition function Z has a remarkably 
simple expression: 

Z= y^6yt with: t= sinh(2 1 J\/kT) (20) 

and where y is the solution of the algebraic equation: 

( ? y-l) 3 (l + 3y) (l + t 2 )-t = 2(l-t) 4 -y 3 (21) 

The partition function has a singular part which behaves as \t — 1| 4 / 3 . 

Some interesting duality properties of the Baxter- Wu model are very clearly 
detailed in p5| , and used to convert the Baxter- Wu model into a coloring prob- 
lem. This provides a very heuristic example showing that duality is not specific 
to edge- interaction spin models, but can also be introduced with many-body 
interactions. In the following we briefly describe how the Baxter- Wu model is 
transformed into a coloring problem. 

First we introduce a Z2-Fourier transform with function g(X, /1) which en- 
ables us to write simply the Kramers- Wannier duality for this three-spin model 
(A and /1 are Ising spins) as: 

g(X, At) = +1 if: A = +1, 

g(X, /x) = fi if: A = -1. (22) 

Note that this function is symmetric in A and /1 , namely, g(X, fi) = +1 when 
At = +1 and <?(A, fi) = X when /i = — 1. 

Returning to the Baxter- Wu model, each spin o~i of the triangular lattice 
belongs to six triangles around vertex i, which form a hexagon with the spin 
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<Zj at the center. Let us now consider the close-packing of such hexagons. The 
spins <7j now form a (triangular) sublattice of the initial triangular lattice. 

Consider next the spins <7j and denote the edge connecting nearest-neighboring 
spins <7fe and 07 sitting on the hexagon surrounding <jj by < kl > . Let us 
introduce Ising edge variables A r corresponding to the six edges < kl > of 
the hexagon: X r = • 07, r = 1, • • - 6. The local Boltzmann weight of a 
hexagonal cell around a spin cr, can be written as: 

6 6 

w hcx = - • (1 + n a -) • ex p(^ • ^ ■ n a ( 23 ) 

r — 1 r — 1 

where the factor (1 + J| A r ) takes into account of the fact that the Ising edge 
variables A r are not independent but constrained by the condition Y[ -V 
1 , As usual this condition, associated with every hexagon, can be written by 
introducing a dummy variable \Xi also associated with every hexagon: 

6 6 

II ^ Mi) = 1 + II K ' 

Hi = ±1 r = 1 r — 1 



enabling us to rewrite (|23|) as: 

6 6 

Whc X = TT 9(k, w) • exp(V ■ a,- ■ A, 

|H = il r = l r — 1 

The partition function of the Baxter- Wu model is now seen as a summation over 
all the (initial) spins <7i and the (dummy) spins /i.; of a triangular sublattice, 
and the edge Ising spins A r . Let us focus on one A r . The edge r — < kl > 
belongs to a hexagon around spin <7j and a neighboring hexagon around another 
spin, say, Oj . The edge Ising spin A r thus occurs in the Boltzmann factors with 
a factor: 

W Xr = e K <°* +^' A " ■ g{\ r , Mi ) • g(\ r , h)- 

Summing over the edge Ising spin A r in the partition function and using relations 
(p2|), one thus obtains a factor: 

= E = e K{pi+ai) + ^^-e- K ^+^ (24) 

x r = ±1 

between two spins tr^ and cr^ on the sublattice. This can be interpreted as the 
edge weight associated with a coloring problem, and the Baxter- Wu model is 
transformed into a coloring problem. 

4.3 The Blume-Emery-Grifflths Model 

The Blume-Emery-Griffiths (B EG) mode l is a mod el that F. Y. Wu and cowork- 



ers considered quite naturally [106| 116 , |l3q , [148|, sin ce it redu ces to an Ising 



model on the honeycomb lattice in a special manifold [106, 11£] 
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The BEG model is defined by the Hamiltonian: 
-0-H = -J- SiSj-K- 53 S? S 2 - Sf - H ■ ^ (25) 

<i,j> <i,j> i i 

where the spins are classical spin-1 spins taking on the values Si = 0, ±1. In 
the high-temperature expansion the nearest-neighbor Boltzmann factor assumes 
the form 

exp( J Si Sj + K Sj S?) = (26) 
1 + (e K sinh J) S t Sj + (e K cosh J - 1) S 2 S 2 . 
It follows then in the subspace K = — ln(cosh J) , one has the simple relation 

exp(J Si Sj + K S 2 S]) = 1 + Si Sj tanh(J). 
and the partition function of the BEG model assumes the simpler form: 

zbeg = n ( x + Si s j tanh j ) n ex p(- A s i + h s i)- 

Si =0, ±1 <i.j> i 

Expanding the products over neighboring pairs, representing each term by a 
graph and making use of the identities 

£ Sf-e" AS ? = p(n) with: p(2)=2 e - A , 

s,= o,±i 

p(0) = 2e- A + 1, p(l) = 0, (27) 

one finds that one has eight possible configurations at each vertex of the honey- 
comb lattice, corresponding to the following Boltzmann weights of an (isotropic) 
eight-vertex model: 

a = 1 + 2 e- A cosh(ff), b = 2 ^tanh(J) e" A sinh(iJ), 

c = 2 tanh(J) e - A cosh(iJ), d = 2 (tanh(J)) 3/2 e~ A sinh(iJ) 

With these notations one deduces the identity of the partition function of the 
BEG mo del (p5|) w ith the eight- vertex model on a honeycomb lattice (see also 
sections (5.1) and ( |5.3|) below): 



J BEG 



Z Sv (a, b, c, d) 



Performing a weak-graph duality transformation on this eight-vertex model (see 
section (||) below) associated with the 2x2 (gauge) matrix Q] : 

i y 



91 = 92= T2 K-v i 

one finds that the partition function of the eight-vertex model Z Sv ( a, 6, c, d) 
remains invariant under a weak-graph duality transformation [Q, 116]: 

~a = (a +3yb +3y 2 c+y 3 d)/(l + y 2 f j2 , • • • 

Z 8v (a, b, c, d) = Z 8v (a, b, c, d) (28) 
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The four parameters a, 6, c, d or a, 6, c, d can be seen, as far as the calculation 
of the partition function is concerned, as four homogeneous parameters. Taking 
into account an irrelevant overall factor and the irrelevant gauge variable y 
from the weak-graph symmetry ( |28| ) , one sees that the partition function of the 
eight-vertex model Z$ v (a, 6, c, d) basically depends on two variables instead 
of four. Not surprisingly, Wu found that Zg v (a, b, c, d) is equivalent to the 
partition function of an Ising model with nearest-neighbor interactions Kj and 
a magnetic field L: 

(L, Kj) = ]~[ exp(i^/ Gi o-j) JJexp(L a t ) 

£T <ij> i 

„ 7 _ 7n /2 cosh(i) cosh 372 ^/)^ , . 
= Z 8 „(a, 6, c, d) ■ [ (29) 

where N is the number of lattice sites. The explicit expressions of Kj and L in 
terms of the BEG parameters are complicated. But for H = 0, one has L = 
and 

2 

tanh(ifj) = r- • tanh(J), 

2 + 

using which one determines the critical line Kj = in the J > regime. 

The spontaneous magnetization of the BEG model for J > and the phase 



boundary of the J < BEG model can be similarly determined [ 1 1 6 1 . 

The proof of the equivalence of the honeycomb eight-vertex model with an 
Ising model in a field as outlined in the above is quite tedious. However, a more 
direct derivation has since been given by Wu ]130| . 

The eight-vertex model on the honeycomb lattice can also be seen to be 
related to a lattice-gas grand partition function "E. Ka9 (z, J, J3) on the Kagome 



lattice [ 126 1 (see also section (6.3) below), where —J is the nearest-neighbor 
interaction, — J3 the triplet interactions existing among three sites surrounding 
a triangular face of the lattice, and z denotes the fugacity. Then one has the 
equivalence: 

Z Ka °(z, J, J 3 ) = Z 8v (a, 6, c, d) (30) 



a 



1, b = V£, c= z-e J , d= z 3 / 2 -e 3J+j3 . 



From (|29|) and (|30|), F. Y. Wu and X. N. Wu were able to obtain results for the 
liquid and vapor densities, showing that an observed anomalous critical behavior 



occurs in the lattice gas only when there are nonzero triplet interactions [ 126 1 



This analysis has been extended to a lattice gas on the 3—12 lattice by J. L. 



Ting, S. C. Lin and F. Y. Wu [140| 



Wu's tricks for the honeycomb BEG model are not limited to the weak- 
graph transformation for the eight-vertex model. Using a syzygy analysis of 
the invariants under the 0(3) transformation L. H. Gwa and F. Y. Wu have 
obtained an expression for the critical vari ety of the hon eyco mb BEG model to 
an extremely high degree of accuracy [145 ] (see section ( |5.4|) below). 
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4.4 Other spin results: disorder points 



Let us finally describe, among many results obtained by F. Y. Wu on spin mod- 
els, one result concerning disorder points. Disorder solutions are particularly 
simple solutions corresponding to some "dimensional reduction" of the model, 
which provides simple exact results for models which are generically quite in- 
volved. While this yields severe constraints on the phase diagrams, the series 
expansion and the analyticity properties of the model, it does lead to exact 
solutions of models which are otherwise nonintegra ble. 

For example, using a decimation approach, Wu [IOC] has deduced the disor- 
der solutio n for the triangular Ising model in a nonzero magnetic field. Wu and 
K. Y. Lin [120 1 have used a checkerboard Ising lattice to illustrate that there 
may exist more than one disorder points in a given spin system. Along the same 
vine N. C. Chao and Wu [ 101 1 have explored the validity of the decimation ap- 
proach by considering disorder solutions of a general checkerboard Ising model 
in a field. 



5 Weak-graph dualities and Hilbert's syzygies 

In a pioneering paper F. Y. Wu and Y. K. Wang |58) introduced a duality 
transformation for a general spin model which can have chiral interactions. This 
is the first time that chiral spin model has been explicitly considered. In terms 
of R matrices such as ([jj]) these transformations are the tensor product of two 
similarities: 

R — ► 9l ®g 2 R ■ g^^g^ 1 (31) 

where g\ and gi are two q x q matrices, R is a q 2 x q 2 matrix ( q — 2 for 
the sixteen- vertex model). This symmetry group is an s\(q) x sl(q) symmetry 
group. The high- and low-temperature duality (||) given in section ( [3.l| ) for the 
sixteen- vertex model is a particular case of such transformations, corresponding 
to <?i and <?2 being two involutions: 

1 / 1 i 
91 = 92= W\~i 1 



This duality relation [p8L 131 is now well-known for vertex models. It corre- 
sponds to symmetries of the model and can be used, as will be seen in the 
next two subsections, to find good variables to express the critical manifolds 
of a lattice model, and, hopefully, to determine their exact expressions when 
algebraic. 



5.1 Hilbert's syzygies, gauge-like dualities and critical man- 
ifolds 

Let us give some hints as to how the gauge-like dualities enable us to deduce 
results on critical manifolds or varieties. The main idea is to construct algebraic 
invariants under these gauge transformations. 
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Hilbert has shown that all invariants of a linear transformation are algebraic 
and can be expressed in terms of a set of homogeneous polynomials, the syzy- 
gies. Considering the sixteen- vertex model, the transformation is 0(2) and the 
fundamental invariants corresponding to the 0( 2) gr oup have been constructed 



by J. H. H. Perk, F. Y. Wu and X. N. Wu in |131|. Likewise for 3-state ver 



tex models the transformation is 0(3) and the associated invariants have been 



constructed by L. H. Gwa and F. Y. Wu [ 146 1 (see section (5.4) below). 



5.2 Hilbert's syzygies and the square lattice Ising model 
in a magnetic field 

With an algebraic prejudice for critical manifolds, it is very tempting to con- 
jecture closed algebraic formula for critical manifolds that will be reproduce 
known exact results in various limits. For instance, closed-form expressions for 
the critical line of the square lattice antifcrromagnctic Ising model in a magnetic 
field were proposed^ by Muller-Hartmann and Zittartz. However, it has been 
shown that the expression is numerically incorrect. 



X. N. Wu and F. Y. Wu [ 135 1 considered the square lattice antiferromagnetic 
Ising model in a magnetic field which can be seen as a subcase of the sixteen- 
vertex model under the 0(2) group. Introducing the variables: 

a = 1, b = \Jvh, c = v, d = v 3 ^ 2 h, e = v 2 
with: v = tanh( J/kT), h = ttmh{H/kT), 

the five fundamental Hilbert's invariants of 0(2) read: 

h = a + 2 c + e, h = (a - 6 c + e) 2 + 16 (b - df , 

h = (a - e) 2 + 4 (b + d) 2 , h = a 2 d - be 2 - 3 {a - e) (b + d) c 

h= (a - 6 c + e) • ((a - e) 2 - 4 (6 + d) 2 ) + 4 (a - e) (b 2 - d 2 ) 

and the critical line proposed by Wu and Wu assumes the form: 

Cl • If C 2 ■ 1\ I 2 + c 3 • If h+ C4 • II + c 5 • If + ce • h h = 0. 

They then determined the c^'s using the various known results including a 
constraint dictated by the known zero- field critical point, as well as the results of 
finite-size analyses which they carried out. This leads to the values c\ = 1, ci = 
-0.044 338, c 3 = 0.362 73, c 4 = 0.000 4938, c 5 = 0.042 779, c 6 = -0.008 9149 
The resulting closed form expression for the critical line reproduces all known 
numerical data to a high degree of accuracy. For instance, the critical line 
yields for small magnetic field: T c ~ T -(l - u-(H/J) 2 ) with u ~ 0.038 022. 
This is compared to the presumably exact value obtained by M. Kaufmann: 
u ~ 0.038 0123 259 •• • 



15 This is just one example in a very long list of wrong algebraic conjectures for critical 
manifolds one can find in the literature. 
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5.3 Hilbert's syzygies and the honeycomb lattice Ising 
model in a magnetic field 

Similar analyses have been carried out for the honeycomb lattice by F. Y. Wu, X. 
N. Wu and H. W. J. Blote ]132fl. For the corresponding honeycomb eight- vertex 



model we have: 

a = 1, b = \/v h, c = v, d = v 3 ^ 2 h 

with: v = tanh(J/fcT), h = t&nh(H/kT). 

Analogous to (|3^), we introduce the following Hilbert's syzygies: 

P=a 2 +3ac+3bd+d 2 i Q = b 2 - a c + c 2 - b d 

P 2 = 2 (a 4 + d 4 ) - 6 (a 2 c 2 + b 2 d 2 ) + 12 (a 2 b 2 + c 2 d 2 ) - 5a 2 d 2 
+ 27 b 2 c 2 + 36 (a 6 + cd) be + 18 abed. 



The critical line proposed by Wu, Wu and Blote now reads 125]: 

Ci-P 2 +c 2 -P 2 +c 3 -PQ + c A -Q 2 = 0. (32) 

After an extensive search by mapping with all known exact results, they pro- 
posed the numbers: c\ = 1, c 2 = —(4 + 3 \/3)/6, C3 = —(1 — 9 \/3)/8, and 
c 4 = -3(3 - V3)/8. 

The initial slope of this critical frontier for small H is — ln(z c ) where z c 
is the critical fugacity of the nearest-neighbor exclusion gas. Their expression 
leads to the value z c ~ 7.851 780 04 • • • which is in very good agreement with 
the value obtained from finite-size analysis, namely, z c ~ 7.851 725 175(13). 
The critical line ([32|) is probably not the exact one but certainly a very accurate 
approximation. 

Comment: The Hilbert's invariant theory amounts to considering linear 
gauge-like symmetries of the model and the associated invariants. From the 
inversion relations one has further an infinite discrete set of birational non- 
linear symmetries, one can couple with these continuous linear groups. In fact 
all the above analyses can be revisited by combining the gauge transformation 
with the infinite discrete symmetries generated by the inversion relations of 
the sixteen- (or simply eight)-vertex models. This would lead us to consider a 
unique "superinvariant" which is, in fact, the modular invariant of the elliptic 
curves parametrizing the sixteen-vertex model. 



5.4 Hilbert's syzygies and the honeycomb BEG model 

To apply the syzygy consideration to the BEG model which is a 3-state spin 
model, one needs to consider the 0(3) gauge transformation and its associated 
invariants. But the construction of the 0(3) invariants is very complicated. 
However, the day is saved since there exists a mapping between 0(3) and sl(2), 
and invariants for the latter have been worked out by mathematicians a long 
time ago. While the mathematics to decipher the old results is involved, L. H. 
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Gwa and F. Y. Wu Il46l have succeeded in carrying out such an analysis and 
deduced that there are 5 independent invariants for the 0(3) group. They nex t 
applied the analysis to the isotropic BEG model on the honeycomb lattice [ 148 1 , 
and found one of the invariants to vanish identically. The remaining 4 invariants 
are then used to determine the critical variety of the BEG model as in the case 
of the 0(2) gauge. The resulting closed-form expression for the critical variety 
agrees extremely well with numbers obtained from a finite-size analysis, which 
they also carried out. 



6 Critical manifolds and critical varieties 

A problem solver like F. Y. Wu first tries to find the exact solution of a problem. 
He tries to "dig out" problems that can be solved. However since most of the 
problems one looks at cannot be solved exactly, one then tries to study models 
for which some exact results can be "salvaged" . It could be the critical manifolds 
which are submanifolds along which the models are Yang-Baxter integrable. 
In such cases the critical manifolds are, in fact, critical varieties. For other 
models the critical manifolds are algebraic varieties without hidden Yang-Baxter 



integrability J79|, |l61 | 



For two-dimensional lattice models the situation is more specific: one can 
have some "conformal prejudice" that critical manifolds should be submanifolds 
where the model has a two-dimensional conformal (infinite) symmetry yielding 
some integrability in the scaling limit. Therefore, as far as critical manifolds 
are concerned, it is crucial to understand the inter-relation between 1) alge- 
braicity consequence of Yang-Baxter integrability, 2) conformal integrability 
consequence of a two dimensional criticality, and 3) self-duality. 

Many criticality conditions have been obtained, or simply conjectured, in 
the literature of lattice models in statistical mechanics jlT], |6^, |72|, [73| |74|, [79], 



|, and all these conjectures were algebraic [ 125 1 . A straightforward situation 
corresponds to the case where the model possesses a duality symmetry (see 
section (^J)) for which it is always possible to give a linear representation of this 
duality transformation. One can sometimes find varieties which are globally 
invariant under this symmetry. 

Let us, instead, consider the fixed points of the linear duality transformation, 
which belong to some algebraic variety (hyperplane) . If the algebraic variety 
separates the phase diagram into two disconnected parts, and if one assumes 
that the critical temperature is unique, one can actually deduce that this alge- 
braic variety is a critical variety |79| . Of course if the algebraic variety is only 
globally invariant (and not invariant point-by-point on the algebraic variety) one 
cannot draw any conclusion. 

In fact, for most of the time one is not in a situation where a simple self- 
dual argument allows the determination of the critical points. A good example 
is the duality transformation ( |3l|) for the sixteen-vertex model for which self- 
dual arguments are insufficient. In fact there exists a "super-invariant" on this 
model after taking into account the gauge (weak-graph) duality symmetries ( |3l| ) 
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and the inversion relation symmetries^. But that is another story. 

When the critical varieties are exact, they are almost always related to some 
integrability of the model, the algebraicity thus being a consequence of the in- 
tegrability. A paradigm is the standard scalar Potts model (see section (7.1)) 
which is integrable at criticality. However, we will give below, in the case of 
the two- and three-site interaction Potts model [|l6l| , an algebraic variety which 



is the critical condition [125] but is unrelated to any simple Yang-Baxter-like 
integrability. 

Throughout the years F. Y. Wu has obtained numerous results of critical 
varieties for two- and three-dimensional spin models, and many related conjec- 
tures as well. One can only say that the seeking of critical manifolds and critical 
varieties is a fascinating subject matter by itself for specialists like Wu. 



6.1 Inversion relations, duality and critical varieties 

The considerations of inversion relations has been shown to be a powerful tool 
for analyzing the phase diagram of lattice models and, particularly, for obtaining 
critical algebraic manifolds in the form of algebraic varieties (see ( |5^ ) in section 
©)• 

Let us consider the standard scalar g-state Potts model on an anisotropic 
triangular lattice with nearest-neighbor and three-spin interactions around up- 
pointing triangles (79 



161] as shown: 




16 Remarkably this consideration even extends the s/(2)®si(2) weak-graph symmetry group 
to a sl(2) (g) sl(2) (gi sl(2) eg) sl(2) symmetry group. 



2cS 



The partition function of the models reads: 



Z = II e Kl J] e K25 ° j '" k I] e*" 3 " 5 ^ T\ e K S ^"i s "i-K 

{a-i} <i,j> <3,k> <k,l> A 

Here the summation is taken over all spin configurations, the first three products 
denote edge Boltzmann weights and the last product is over all up-pointing 
triangles. 

A duality transformation exists for this model . We introduce the follow- 
ing notation: 

— & , — c ' . i — _L . 2.- . • > 

= x Xi x 2 x 3 — (xi + x 2 + X3) + 2 (33) 




q 2 



y 

x x + x 2 + x 3 - 2 + q 2 /y 



(34) 



1 ^2 ^3 



and the partition transforms as 

Z(x!,x 2 ,x 3 ,y) = {y/q) N ■ Z{x\,x%,xl,y*) (35) 

where N is the number of sites. 

On the basis of this duality Baxter et al. proposed that the critical points 
are located on the algebraic variety: 

x xi x 2 x 3 - (xi + x 2 + x 3 ) + 2 - q = (36) 

which corresponds to the set of fixed points of D. The critical variety (^) is 
not only globally invariant under (|34|), it is also point-by-point invariant, namely, 
every point on the variety is invariant. In general when an algebraic variety is 
such that every point of the variety is invariant under a duality symmetry, it is 
possible to argue, subject to some continuity and uniqueness arguments, that 
the variety actually corresponds to the criticality variety. This has been done 



by Wu and Zia 125] for q > 4 in the ferromagnetic region. It is important to 



note that the critical variety (36) is not an algebraic variety on which the model 
becomes Yang-Baxter (star-triangle) integrable. This is an interesting example 
of a model where algebraic criticality does not automatically imply Yang-Baxter 
integrability. 

Comment: In suitable variables the duality transformations can be seen as 
a linear transformation. There are two globally invariant hyperplanes under D: 
y = +q and y = —q. The (ferromagnetic) criticality variety ( |36"1 ) corresponds 
to y = +q. The second hyperplane y = —q is not a point-by-point invariant 
although it is globally self-dual. It is not a locus for critical or transition points. 
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This illustrates a fundamental question one frequently encounters when try- 
ing to analyze a lattice model: is the critical manifold an algebraic variety or a 
transcendental manifold? It will be seen that a first-order transition manifold 
exists for this model for q = 3, and its algebraic o r tran scendental status is 
far from being clear (see |l66| and j67| ) in section ( 7 . 4 . 3| ) ) . The existence of 



such very large (nonlinear) group of (birational) symmetries provides drastic 
constraints on the critical manifold and therefore the phase diagram. 

There exist three inversion relations associated with the three directions of 



the triangular lattice exist for this model [161|. For instance, the inversion 
relation which singles out direction 1 (see figure 1) is the (involutive) rational 
transformation I\. 

( (mi - l) 2 (zi +q-2) 

h ■ [X, Xl, X2, xs) — • 



V (x x\ + x X\ (q — 3) — q + 2 ) (x\ — 1) ' 

xi [x - 1) xi - 1 x x - 1 \ , , 

2-q-X! + — , — — , — — ). 37) 

X\X—\ X-i\XX\~ 1) X2\XX\ — 1) / 

These three inversion relations generate a group of symmetries which is naturally 
represented in terms of birational transformations in a four dimensional space. 
This infinite discrete group of birational symmetries is generically a very large 
one (as large as a free group). The algebraic variety (f56|) is remarkable from 
an algebraic geometry viewpoint: it is invariant under this very large group 
generated by three involutions (|37j). 

In this framework of a very large group of symmetries of the model, an 
amazing situation arises: the one for which q, the number of states of the Potts 
model, corresponds to Tutte-Beraha numbers q — 2+2 cos(27r / TV) where iV 
is an integer. For these selected numbers of q the group of birational transfor- 
mations is generated by generators of finite order: it is seen as a Coxeter group 
generated by generators and relations between the generators. The elements of 
the group can be seen as the words one can build from an alphabet of three let- 
ters A, B and C with the constraints A N+1 = A, B N+1 = B, C N+1 = C. 
Since the generators A, B and C do not commute (and nor do any power of 
A, B and C) the number of words of length L still grows exponentially with L 
(hyperbolic group). Among these values of q, two Tutte-Beraha numbers play 
a special role: q — 1 and q = 3. For these two values the hyperbolic Coxeter 
group degenerates^] into a group isomorphic to Z x Z. 

For the standard scalar nearest-neighbor Potts model the Tutte-Beraha num- 
bers correspond to the values of q for which the critical exponents of the model 



are rational (see (p3|) in section (7.1)) 



Up to semi-direct products by finite groups. 
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6.2 The exact critical frontier of the Potts model on the 
3-12 lattice 

F. Y. Wu et al. considered a general 3—12 lattice with two and three-site inter- 



action on the triangular cells [155]. This model has eleven coupling constants 
and includes the Kagome lattice as a special case. 

In a special parameter subspace of the model, condition (|3^) below, an exact 
critical frontier for this Potts model on a general 3—12 lattice Potts model was 
determined. The Kagome lattice limit is unfortunately not compatible with the 



required condition (|38| ) 



The condition under which they obtained the exact critical frontier reads: 

x 2 x\ x\x\ — x xi X2 x 3 ■ (xi X2 + x 2 x 3 + xi x 3 — 1) (38) 
+ {xi + x 2 + X3 + q - 4) • (xi x 2 + x 2 x 3 + x\ X3 + 3 — q) 
— qxi x 2 x 3 — {x\ + x\ + £3) + q 2 — 6 q + 10 = 

This is nothing but the condition corresponding to the star-triangle relation of 
the Potts model. 



Comment: One can show that condition ( |38] ) is actually invariant under the 
inversion relation ( |37| ) of the previous section (3.1), and therefore since (38) is 



symmetric under the permutations of K\, K 2 and K 3 under the three inver- 
sions generating the very large group of birational transformations previously 
mentioned in section (|6.l[ ). More generally, introducing Di, D 2 and D 3 : 

Dx — xi + x 2 + x 3 - x x\ x 2 x 3 + q - 2 , D 3 — x X\ x 2 x 3 - xi x 2 x 3 
D 2 = x\ + x 2 + x 3 + xxix 2 x 3 — 1 — (xi x 2 + x 2 x 3 + xi x 3 ), 

one can show that the algebraic expression 

Ii(xi, x 2 , x 3 , x) = — — ^ D2 (39) 
L>i L> 2 — q ■ L> 3 

is invariant under the three inversion relations and the large group of birational 
transformation they generate, the (star-triangle) condition fl38| ) corresponding 
to l\(xi, x 2 , x 3 , x) — 00, namely^] D\D 2 — q D 3 = 0. 

When x = 1, or q = 1 or 3, there are additional invariants of the three 
inversions (|37]). For instance, for x = 1, one can build an invariant from a 
covariant we give below (see (|57|)). For q — 3, introducing 

n /2222 2 2 2 2 22,2,2,2 i\ 

— X\ X 2 X 3 ' I X-^X 2 X 3 X X 2 X X-^ X 3 X X-^ ^2"^3"^ ' ^1 ' ^2 ' ^3 / 

one finds that the expression: 

Df-D 2 



I 2 {xi, x 2 , x 3 , x) 



3 5 • D 5 



18 For x = 1 (no three-spin interaction, Dg = 0), condition ( |38| ) factorizes and one recovers 
the ferromagnetic critical condition (^) of the g-state Potts model on an anisotropic triangular 
lattice. 
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is invariant under the three inversions (|37|). One can try to find the mani- 
fold corresponding to the first order transition (see ( 7.4.3| ) below) in the form: 



F(l±, I2) = 0. It still remains an open question to determine whether this va- 
riety is algebraic or transcendental. The x — 1 limit corresponds to Ii = +1. 
The condition 12(2:1, X2, X3, x) = 1 yields X\ — X2 — £3 = 0.215 816 (to 
be compared with 0.226 681 from ( |57|) in section (7J2) below) but still differ- 



ent from 0.204 (see (p6J) in section (7.4.3) below) which is believed to be the 



location of the first-order transition point. 

6.3 The embarrassing Kagome critical manifold 

At the end of the 80's there was a surge of interest in the Kagome lattice com- 
ing from the theoretical study of high-T c or strongly interacting fcrmions in two 
dimensions (the 2D Hubbard model, resonating valence bond (RVB), ground 
state of the Heisenberg model). But the two-dimensional Gutzwiller product 
RVB ansatz strategy promoted by P. W. Anderson for describing strongly in- 
teracting fermions seemed to fail for regular lattices (square, triangular, ...). 
Thus, because of its ground state entropy and other specific properties, the 
Kagome lattice seemed to be the "last chance" for the RVB approach. 

Since one can obtain a critical frontier ( |38| ) for the general 3—12 lattice 
model, and since the 3—12 model includes the Kagome lattice as a special case, 
it is tempting to try to obtain the critical frontier for the Potts Kagome lattice. 

The Kagome Potts critical point was first conjectured by Wu [Q as: 

y 6 - 6 y 4 + 2 (2 - q) ■ y 3 + 3 (3 - 2 q) ■ y 2 - 6 (q - 1) • (q - 2) y 

-(q-2) (q 2 -4 q + 2) = (40) 

which gives for q = 2 the correct critical point y 4 — 6 y 2 — 3 = and for 
q = gives (also correctly) y = 1. Furthermore, for large q, y behaves like 
yjq as it should. However in the percolation limit q — > 1, it gives a percolation 
threshold p c for the Kagome lattice of: p c = 0.524 43 • • • which is compared to 
the best numerical estimate^ obtained by R.M. Ziff and P.N. Suding, namely 
p c = 0.524 405 3 • • • with uncertainty in the last quoted digit. Wu's conjecture 
is thus wrong but by less than 5 • 10~ 5 . Some very long high-temperature series 
of I. Jensen, A. J. Guttmann and I. G. Enting on the q-state Potts model on 
the Kagome lattice further confirm that the conjecture is wrong for arbitrary 
values of q. Nevertheless the Wu conjecture remains to be an extraordinary 
approximation. 

It is a bit surprising that no exact result on integrability (along some al- 
gebraic subvariety) or exact expression for the critical variety is known for the 
standard scalar Potts model on the Kagome lattice, as generally one expects the 
integrability on one lattice, say the square lattice, implies integrability for most 
of the other Euclidian lattices. This is certainly not the case for the Kagome 
lattice. 



19 R.M. Ziff and P.N. Suding, Determination of the bond perc olation threshold for the 
Kagome lattice, J. Phys. A 30, 5351 (1997) and |cond-mat/9707110| . 
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7 Potts models 



The Potts model encompasses a very large number of problems in statistical 
physics and lattice statistics. The Potts model, which is a generalization of 
the two-component Ising model to q components for arbitrary q, has been the 
subject matter of intense interest in many fields ranging from condensed matter 
to high-energy physics. It is also related to coloring problems in graph theory. 

However, exact results on the Potts model have proven to be extremely 
elusive. Rigorous results are limited, and include essentially only a closed-form 
evaluation of its free energy for q = 2, the Ising model, and critical properties 
for the square, triangular and honeycomb lattices |7(J. Much less is known 
about its correlation functions. 



7.1 Wu's review on the Potts model 

F. Y. Wu's 1982 review on the Potts model is very well-known |89| (see also Q). 
It is an exhaustive expository review of most of the results known on the Potts 
model up to 1981, a time when interests on the model began to mount. It has 
remained extremely valuable for anyone wishing to work on the standard scalar 
Potts model. In particular, it explains the q — > 1 limit of the percolation 
problem (see also the q — > 1/2 limit of the dilute spin glass problem, 

and the q — > limit of the resistor network problem, the equivalences with 



the Whitney- Tutte polynomial p9] (see section (7.7) and also pTJ)) and many 



other related models are also detailed. For instance, the Blumc-Capel and the 



Blume-Emery-Griffiths model (see (£5j) in (4.3)) can also be seen as a Potts 
model. More generally it is shown that any system of classical q-state spins, 
the Potts model included, can be formulated as a spin (q — l)/2 system. 

However, Wu's review was not written in time to include discussions of the 
inversion functional relations. For the two- and three-dimensional anisotropic 
g-state Potts models, the partition functions satisfy, respectively, the functional 
relations: 

Z(e K \ e K2 ) ■ Z(2 - q - e K \ e~ K *) = (e Kl - 1) • (1 - q - e Kl ), (41) 



Z CU Mc(e K \ e K \ e K *) ■ Z cublc (2 - q - e K \ e~ K \ e Ks ) 

= {e Kl -l)-(l-q -e Kl ). (42) 

There are also permutation symmetries like, in 3 dimensions, Z cu bi c (e 1 , e K2 , 
e K3 )= Z cuUc {e K \e K \e K2 )= Z cuUc {e K \ e K \ e Kl ) . Combining these rela- 
tions one generates an infinite set of discrete symmetries which yield a canoni- 
cal rational parametrization of the Potts model at and beyondp^] T — T c , and 
shows clearly the role played by the Tutte-Beraha numbers. These infinite sets 
of discrete symmetries impose very severe constraints on the critical manifolds 

20 At T = T c one recovers the well-known rational parametrization of the model (occurring 
in ( p^ ) see below). 



33 



and the integrability (see sections (fTl|)). The inv ersion relation study has 
subsequently been carried out by F. Y. Wu et al. [ 161 1 . 

Graph theory plays a central role in Wu's work on the Potts model. The 
Potts partition function can be written as [ |39"| : 

Z = Z G (q,K) = £ (e K -l) b q n , (43) 

G'CG 

where K — J/kT, the summation is taken over all subgraphs G' C G, and 
b and n are, respectively, the number of edges and clusters, including isolated 
vertices, of G' . The duality relation of the Potts model can be obtained from a 
graph-theoretical viewpoint by using the Euler relation c + N = b + n where 
c is the number of independent circuits in the subgraph G", and N is the total 
number of vertices in G. This leads to the duality relation 

Z G (q,K) = ■ q 1 -^ ■ Z D (q,K*), 

where D is the graph dual to G, and the dual variable K* is given by: 

(e K - 1) • (e K * - 1) = q. (44) 

The generalization of the duality to multisite interactions is also given in the 
review. 

A consequence of (^) is that one finds the following connection with the 
chromatic polynomial Pq (q) on G by taking the antiferromagnetic zero-temper- 
ature limit K — > — oo: 

Z G (q, K = -OO) = P G (q). 

The (high- and low-temperature) series expansions are described from a 
graph-theoretical viewpoint. For instance, the high-temperature expansions are 
written in the (Domb) form: 

9-1 

z G ( q , k)= J2 II — • C 1 + &)> fa = -^rz ■ (<*> - x ) ( 45 ) 

The introduction of these variables comes from the fact that: 

9-1 

E fa = ( 46 ) 

dj =0 

and, consequently, all subgraphs with vertices of degree 1 give rise to zero con- 
tributions. The number of subgraphs that occur in the expansion is therefore 
greatly reduced. 

The location of the critical points of the anisotropic Potts model on a square, 
triangular and honeycomb lattice were given in terms of the variables x r — 
(e Kr — (see also ( |6l"|) below). These expressions are invariably the various 



34 



special cases of Wu's conjecture 89 ) for the critical point of the more general 
checkerboard lattice, namely, 



y/q + x\ + X2 + x 3 +X4 (47) 

= X\ X 2 X 3 + X\ X2 X4 + X\ X 3 X4 + X 2 X 3 X4 + yfq ■ X\ X 2 X 3 X4- 

Using the notation a = e Kl , b = e K2 , c = e K3 , d = e Ki , this critical 
algebraic variety (^) reads: 

-{q-l)(q-3) + (a + b + c + d)(2-q) (48) 
— (ab + ac + bc + ad + bd + cd) + abed = 

The critical point of a mixed ferromagnetic-antiferromagnetic square Potts model 
considered by Kinzel, Selke and Wu |82| : 

(e Kl - 1) • (e K * + 1) = -q (49) 

was also given. While this expression coincides with the critical point for q = 2, 
it is incompatible with the inversion relations ( f4l|) for general q, and hence is 
not a critical variety. This is because the infinite discrete group of symmetries 
generated from the inversion relations of the square Potts model transforms ( ^9|) 
into an infinite set of other algebraic varieties, and hence cannot be critical. 

Generally, critical manifolds need to be (globally) invariant under this infi- 
nite set of transformations (discrete symmetries). Actually, for the anisotropic 
square Potts model, for instance, one can show that, when q is not a Tuttc- 
Beraha number 21 , the only algebraic varieties compatible with the inversion re- 
lation symmetries ( |4l| ) are given by the well-known ferromagnetic condition: 

(e Kl - 1) • (e K2 - 1) = q (50) 

and the antiferromagnetic condition obtained by R. J. Baxter: 

{e Kl + 1) ■ {e I<2 + 1) = 4 - q (51) 

(for which the model is exactly soluble). Note that these two varieties can be 
deduced from the conjecture (E?l) by taking xx = a; 3 and X2 = X4. In this limit 



the critical condition (47) factors into conditions (BOJ) and (plj). In fact, it has 



since been shown that the critical condition (47) corresponds to an integrability 
condition of the checkerboard Potts model. 

At criticality the Potts model is exactly solvable. Let us give the exam- 
ple of the square lattice. The free energy of the isotropic Potts model at the 
ferromagnetic critical point T = T c reads: 

/(*!■.). 1 + , + 2- ± <- • •'"""■"> , fc: 9 >4 

2 * — ' n 

71=1 



21 When q is a Tutte-Beraha number the generically infinite discrete group T generated 
by the inversion relations on the model becomes a finite group and many algebraic varieties 
can be invariant under such a finite group: a simple way to build such algebraic T-invariants 
amounts to performing summations over the group T. 
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f(q, T c ) = ln(2) + 4'ln(^|) for: q = 4 (52) 

r, m \ 1 dx , . . sinh((7T — it) x) 

/(g, T c ) = - + / — tanh (/x as) 7 for: 9 < 4 

2 ,/_oo a; smh(7r x) 

where the variables and correspond to the rational parametrization of the 
model at T — T c , namely cosh(#) = or cos(fi) — y/q/2. 

The exact critical exponents of the standard scalar Potts model are also 
given in this review. These critical exponents, which are rational when q is the 
Tuttc-Bcraha numbers, are: 

, 2 1 - 2 u . 1 + u , 7 - Au + u 2 

a = a = - ■ - , (3 = — — , 7 = 7 



3 1 - u ' ' 12 ' ' ' 6-(l-ti) 

(3 - u) • (5 - it) , 2-u 1-u 2 
S = 1 5 > v = v = — — -, r\ = 



1-u 2 ' 3 -(1-u)' ' 2 -(2-u) 

where the parameter u is related to q by: 

J = \fq, or: 2 + 2 cos(7ru) = q (53) 

These results play a key role in the emergence of the conformal theory. 

7.2 Comments on the checkerboard Potts model 

The Wu conjecture (|47f) for the criticality condition of the g-state checkerboard 
has since been confirmed from an inversion relation analysis. To discuss the 
inversion relation we introduce variables u, v, w, z, t defined by: 

e Kl - 1 _ u-t e K2 - I _ v-t 

e A 'i +q-l ~ *' l-t 3 u W* +q- 1 ~ * ' 1 - t 3 V ' "' 11 

or: = a - .*-».- U - *" 



t- (1 -it*)' i< (1 

with: i + = y/q 

In these variables the criticality condition ( ft7| ) reads: 

uvwz = 1. (55) 

Using the inversion trickQ the partition function of the checkerboard model at 
criticality can be written in a multiplicative form 

Z(u, v, w, z) [uvwz = i] (56) 
g F(it) -F(l/u) F(w)-F(l/w) F(w)-F(l/w) F(z)-F(l/z) 
t 2 1 - tu 1 - tv 1 - tw 1 - tz 



where: F( u ) = j [ j 



1 - t 



4 n-1 , 



f4 n+1 , 



22 With the well-suited variables ( |54[ ) the inversion trick naturally yields exact formulae as 
products over an infinite discrete group, in this case Eulerian products (p6[) - 
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This formula is reminiscent of ( |Io| ) for the six- vertex model (see section ( |3.l| )), as 
it should since it also follows from the Potts and six- vertex model correspondence 



(see section ( |7.3| ) below). 

Note that, in the anisotropic case w — u and z = v, the critical condition 
( p5|) factorizes into mo = +1 and nti = -1 which are respectively the afore- 
mentioned ferromagnetic and antiferromagnetic criticalitics (and integrability). 

The checkerboard model reduces to an honeycomb model, and hence a trian- 
gular model by taking the dual, if one of the four interactions vanishes. There- 
fore it is useful to re-examine the triangular lattice limit of the checkerboard 
variety. P. Martin et al. have established two varieties for the triangular Potts 
model: a ferromagnetic variety which is precisely (|3^) with x = 1, and an 
antiferromagnetic variety 

(q - 2) 2 - 2 + (a + b + c + abc) (q - 2) + 2 (ac + ab + be) = 0. (57) 

These two algebraic varieties can be written, respectively, as u v w = +t and 
uvw = —t, which can be deduced by taking, respectively, d — 1 and d = 
-(?-2)/2in@. 

For q = 3 the antiferromagnetic algebraic variety yields an algebraic criti- 
cal point v ery close to the first-order transition point for the three-state Potts 
model @,[L43]. 



Along the well-known ferromagnetic variety (|44| ) which for the isotropic 
model reads (e K — 1) • (e K — 1) = q, it is interest to point out the existence 
of some "hidden" duality K — * such that the antiferromagnetic condition 
( |5l| ) reads K2 — K\. This hidden duality is the involution: 

" < 58 > 

Note that the dualities (Q) and ( |58| ) commute, and their product gives the 
involution: 



e K — 



K *t (q-2)-e K +2 

= " 2e*+{ q -2) (59) 



The ferromagnetic critical variety of the anisotropic triangular Potts model ( pq ) 
now transforms into the anti- ferromagnetic critical variety ( |57| ) under (p9|). One 
also finds that the ferromagnetic critical variety of the anisotropic square-lattice 
Potts model (pQ), and its anti- ferromagnetic critical variety (pl|), are both in- 



variant under transformations (44), (pq ) and thus (|5^). More generally, for the 
checkerboard model one finds that ( ff8| ) is invariant under both (^4|), ( pq ) and 
thus (p9|). These results can be very simply seen in the variables u, v, w, z. For 
example, ( |59| ) is simply u — * —u. Involution like the "hidden duality" (pq), or 
like (p9[), do not yield simple functional equations on the partition function like 
the Kramers- Wannier duality does (see (jij) and (|3|)). They are not symmetries 
of the model, but, rather, "symmetries of the second kind": symmetries of the 
symmetries. 
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7.3 Equivalence of the g-state Potts model with a six- 
vertex model 



Recall that the partition function of the standard scalar q-state Potts model 
can readily be written as a Whitney- Tutte polynomial. In a further step Tem- 
perley and Lieb used operator methods to show that, for the square lattice, the 
Whitney- Tutte polynomial is, in turn, equivalent to a staggered ice-type vertex 
model. 

In a classic paper (67) by R. J. Baxter, S. B. Kclland and F. Y. Wu, this 



equivalence is rederived from a graphical approach (see (|72| ) in section (7.7)), 
which is easier than the algebraic method of Temperley and Lieb and is appli- 
cable to arbitrary planar graph. Particular, the equivalence was extended to 
triangular or honeycomb Potts models and a staggered six-vertex model on the 
Kagome lattice. 

The equivalences are as follows (page 404 in J57|). For the square, triangular 
and honeycomb lattices, the equivalent ice- type vertex model has the Boltzmann 
weights: 

{u>i, 0J2, we) = (1, 1, x T , x r , A r , B r ) (60) 

where {A r , B r ) are respectively, for the square, triangular and honeycomb lat- 
tice: 

x 



/ 1 X r \ / 1 , X r \ ( 1 

(-+x r .s,- + s), (- +Xr .t, ¥ + t), ^ 



with: 



s = e 



t = e e / 3 , 2 cosh(6») = Jq, x r = — (61) 

v 7 ? 



These mappings play an important role in the Lee- Yang theorem to be discussed 



latter (see section (7.5.1)) 



7.4 Miscellaneous Potts model phase diagrams 
7.4.1 Potts model with competing interactions 

Kinzel, Selke and Wu |83] have studied a square lattice Potts model with com- 



peting interactions alluded to earlier in section (7T). A similar model in 3 
dimensions with next-nearest-neighbor competing interactions has been studied 
by J. R. Banavar and Wu using the mean-field theory and Monte-Carlo simula- 
tions |)7j . A rich phase diagram was found and they established positively that 
the behavior of the 4-state three dimensional Potts model is mean-field-like. 

7.4.2 First-order transition in the antiferromagnetic Potts model 

F. Y. Wu et al. have analyzed specifically the three-state triangular Potts model 



and considered the (tricky) tricritical behavior of this model |166 | 
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Recalling the results of section ( p.l| ) , and in particular the special role played 
by q = 3 for the two- and three-site interaction Potts model (|33"|), Monte- Carl o 



calculations of the q = 3 isotropic limit of the model have been performed [ 144 1 . 
These studies confirmed the existence of an antiferromagnetic critical point (in 
addition to the well-known ferromagnetic one) probably corresponding to a first- 
order transition occurring near the variety x = 1 isotropic region. 

For the triangular Potts model with two-spin interactions, Monte-Carlo cal- 



culations [ 166 1 confirm the localization a ~ 0.204 ± .003 for a first- order tran- 
sition point. This transition point is confirmed to be different from the algebraic 
antiferromagnetic point localized at a = 0.22665 derived below (see ( |67| ) be- 
low). 

The antiferromagnetic transition point a ~ 0.204 ± .003 can be interpreted 
as to belong to some singular manifold in the parameter space of the model with 
two- and three-spin interactions This singular manifold corresponds to a first- 



order transition frontier. Recalling discussions in section (6.1), the question of 
the algebraic or transcendental status of this first-order transition frontier still 
remains open. 

7.4.3 Chiral Potts models 

F. Y. Wu et al. have analyzed a particular two-dimensional chiral Potts model, 
namely, the 3-state chiral Potts in order to understand the relation between the 



(higher genus) integrability and criticality conditions [142]. On the checkerboard 
model the higher genus integrability of the 3-state chiral Potts model is restricted 
to the following algebraic variety: 

3 • (Qi P 2 P 3 Pi + Q2 Pi P 3 Pi + Q 3 Pi P2 Pi + Qi Pi P2 P 3 ) (62) 
-(Pi Q 2 Q 3 Qi + P2 Qi Q 3 Qi + P 3 Qi Q2 Qi + Pi Qi Q2 Q 3 ) = 

where: 

Pi = fi - 3h h Q t = ft - 2 9i + 3 hi, hi = of bf c\ (63) 
fi = a t h Cl ■ (al + b\ + cf ), 9l = a 3 b\ + b\ c\ + cf a 3 , 

where at, bi, Ci denote the three possible values of the four edge Boltzmann 
weights Wi(o-k — 01) of the checkerboard lattice: 

ai = Wi(Q), h = Wi(l) = Wi(—2), a = Wi(2) = w,:(-l), i = 1, 2, 3, 4 

Recalling the conformal theory prejudice (criticality in two dimensions versus 
integrability), one can also wonder if an integrability condition like (^2|) could 
correspond to a critical subvariety of the phase diagram. Let us consider the 
standard scalar limit of this model. The higher genus integrability condition 
( |62|) reduces to the critical condition ( fi7|) or ( ff8l) of the standard scalar Potts 
model on the checkerboard lattice for q = 3: 

abed - (ab + ad + ac + bc + bd + cd) - (a + b + c + d) = (64) 
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together with another algebraic variety: 



abed + 2 • (acd + bed + abd + abc) + ab + bd + cd + ac + ad + be 

-{a + b + c + d) -2 = (65) 



With the variables (54) taken for q = 3, the critical condition ( |64| ) reads 
u vw z = +1 and the algebraic condition ( |65[ ) reads: uvwz = — 1. Consider- 
ing the similarity of ([35]) , namely uvwz — — 1 , with (|64|) , namely uvwz = +1, 
it is tempting to imagine that (J65j) could also be in some domain of the param- 
eter space a, 6, c, d, a critical variety. 

In the isotropic triangular and standard scalar limit the higher genus in- 
tegrability condition ([32]) factorizes into two conditions. One is the ferro- 
magnetic critical condition of the standard scalar Potts model (see also ([36])) 
2 — q — a — b — c + abc = (or uvw — t) with q = 3, the other one is: 

(q - 2) 2 - 2 + (a + b + c + abc) (q-2) + 2 (ac + ab + be) = 0, 

or: uvw = —t (66) 



with q = 3. In the isotropic limit and for q 
give respectively: 



3, these two algebraic varieties 



3a - 1 







and: 



6 a 2 + 3 a - 1 







(67) 



yielding the ferromagnetic critical point a — 1.8793 and an antiferromag- 
netic transition point at a = 0.22665. This antiferromagnetic point must be 
compared with the antiferromagnetic critical and first-order transition point 
a ~ 0.204 ± 0.003 obtained from series analysis by I. G. Enting and F. Y. 
Wu M. 



7.5 Zeros of partition functions of Potts models 

In 1952 Yang and Lee introduced the concept of considering the zeros of the 
grand partition function of statistical mechanical systems, a consideration that 
has since opened new avenues to the study of phase transitions. While Yang 
and Lee considered the zeros in the complex fugacity plane, or equivalently the 
complex magnetic field plane in the case of spin systems, Fisher in 1964 called 
attention to the relevance of the zeros of the canonical partition function in 
the complex temperature plane. Generally speaking, there exist several differ- 
ent kinds of exact results on lattice models in statistical mechanics. Ideally, 
one would like to obtain the exact, closed- form, expressions of thermodynamic 
quantities such as the per-site free energy, the surface tension, spontaneous mag- 
netization, and correlation functions. A knowledge of these exact expressions 
leads to a complete description of the system including the phase boundary 
(critical frontier) and the location of zeros of the partition function. However, 
exact evaluations of physical quantities are not always possible. In such cases 
one can sometimes determine the critical frontier from properties such as the 
duality and the inversion relations, or analyzes analyticity properties of the free 
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energy by locating the zeros of the partition function. But other than in the 
case of some special one-dimensional model, exact results on the zeros have been 
confined mostly to the Ising model. 



7.5.1 Fugacity variable for checkerboard Potts model, staggering 
field, Lee- Yang theorem, duality 

Let us consider a q-state Potts model on a square, or triangular, lattice, with 
no magnetic field. Since one does not have a magnetic field one does not expect 
a Lee- Yang theorem to exist. Actually this is not true: there is a "hidden" field 
for the standard scalar q-state Potts model ! 

A. Hinterman, H. Kunz and F. Y. Wu |Q combined the equivalence of the 
Potts model with a staggered six-vertex model, together with Lee- Yang circle 
theorem due to Suzuki and Fisher (see also |30|) to deduce the critical variety 
of the Potts model. In this consideration a fugacity variable z is associated 
with the staggering field that occurs in the aforementioned correspondence. In 



terms of the variables {A r , B r ) given in (7.3), the fugacity variable z reads 



respectively for the square and triangular lattices: 

- 4 = - 6 = i^jj 3 - (68) 

B\ D2 B\ B2 B3 

A condition on the Suzuki-Fisher extension of the Lee- Yang theorem and for 
real temperatures requires that we must have q > 4. 

All these results can be generalized to the g-state Potts model on the checker- 
board lattice. Not surprisingly, ( |68| ) is generalized into: 

8 A x A 2 A 3 A 4 

z = 

B\ Bi -B3 -B4 



Note that the duality (44) of the Potts model has a very simple representation 
in terms of these fugacity variables z: z — ► 1/z. 

The zeros of the partition function of the checkerboard model will later be 
seen to lie on \z\ = 1, as it should from the Lee- Yang theorem which applies 
when q > 4. Remarkably, it was also seen later by other authors, that the 
\z\ = 1 condition can be extended to q < 4 and that one recovers the well- 
known Fisher's circles for the Ising model !! 

The fugacity variable z is a fundamental variable. It corresponds to a crucial 
combination of variables encapsulating the action of the infinite discrete group 
generated by the inversion relations. The criticality condition correspond to 
z = 1 and z = —1 seem also to play some role (see (|65|) for q = 3 in 



section (7.4.5)), but not a critical or transition point role. In terms of variables 



u, v, W, z of section ( 7.4.3 ), the fugacity variable z is simply the product uvw z. 



7.5.2 Zeros for the square lattice: A graph-theoretical viewpoint 

Following F. Y. Wu's approach let us consider the g-state Potts model on the 
square lattice from a graph-theoretical viewpoint with the partition function 
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([43|). Introducing the variable 

x = (e K - l)/y/q, (69) 
the partition function ( |4^ ) can be written as a polynomial in x as: 

E 

Z = P G (q,x) = ^c 6 (<z)x & , where c b (q) = q b ' 2 £ <f, 
6=0 G'CG 

where the second summation is taken over all G'CG for a fixed b. Then, the 
duality relation (Q) can be rewritten as a duality relation for the polynomial 
Pa @ 

P G («r,a;) = ^- 1 - E / 2 x E -P D ( g , a; - 1 ). (70) 

In the case of the square lattice for which D is identical to G in the thermody- 
namic limit regardless of boundary conditions, (fro|) implies that the system is 
critical at i c = 1. For fin ite s elf-dual lattices relation ( |70| ) gives an example 
of a self-dual polynomial!^] [ 197 1 . 



To describe the density of zeros on the Lee- Yang circle, we introduce an angle 
9 associated with the location of the zero on the unit circle. For small 9 we have 
g{9) = a^l 1 -"^) f or q < 4 and g{9) = e(q) when q > 4. This leads to the 
specific heat singularity |i|~ Q ( g ) for q < 4 and a jump discontinuity of e(q) in 
U for q > 4. This is the known critical behavior of the Potts model |8i| . 

The zeros of the partition function of the g-state Potts model on the square 



lattice have been evaluated numerically [170 1 . On the basis of these numerical 



results, it was conjectured [ 1 70 1 that, for both finite self-dual lattices and for 
lattices with free or periodic boundary conditions in the thermodynamic limit, 
the zeros in the Re(x) > region of the complex x plane are located on the 
unit circle \x\ = 1. 

7.5.3 Zeros for two- and three-spin interactions on triangular lattice 

We now return to the g-state Potts model on the triangular lattice with two- 
and three-site interactions in alternate triangular faces |n| (section ( |6.l| )). The 
partition function is: 

Z{x,xi,X2,x^) = W(G) , where: 
G 

W(G) = U(l + vS abc )(l + v 1 S bc ){l + v 2 5 ca )(l+v 3 S ab ) 

A 

and v = e K — 1, v,- = e Ki — 1, 

and the product is taken over all up-pointing triangles. It is convenient to 
represent terms in the expansion of the partition function by graphs G in which 

23 A polynomial P(x) in x is self-dual if it is proportional to P(l/x). Self-dual polynomials 
occur naturally in lattice model s in statistical physics and in restricted partitions of an integer 
in number theory (see section (p.3|) below). 
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the up-pointing triangular faces are either occupied by a solid triangle with a 
fugacity v, or unoccupied. 

We next evaluate the weight W(G) associated with the graph G. It is clear 
that each solid triangle contributes a factor v to W(G), and each bond a factor 
Vi. In addition, by including the associated bond factors, each solid triangle 
contributes an additional factor (l + i>i)(l + i>2)(l+V3) = e Kl+K2+Ks . Consider 
next the q dependence of W(G). For the graph representing N isolated points 
we have simply W(G) = q N . For other graphs, each triangle reduces the factor 
q N by q 2 , and each bond by q. But whenever the triangles and bonds close up 
to form a circuit,^] this restores a factor q due to the overlapping of one lattice 
site summation. Thus we have: 



Z = 



E 



,K-l+K 2 +K 3 



- rn 




bi 




b 2 














q 




q 




q 



(71) 



where the summation is over the 2^ graphs G, m is the number of solid tri- 
angles, bi the number of bonds with weight Vi, and c(G) the number of in- 
dependent circuits in G. The expression (Q) generates the high-temperature 
expansion of the partition function. 

In the case of pure three-site interactions, (|7l|) reduces to 



Z 



G 



where 



w 



1)/?- 



For pure three-site interactions, the partition function is self-dual and the critical 
variety assumes the simple form w = 1. On the basis of a reciprocal symmetry 
and numerical results, Wu et al. [170 1 conjectured that zeros for the three-site 
Potts model (in up-pointing triangles) lie in the thermodynamic limit on the 
unit circle \w\ = 1 , as well as a line segment on the real negative axis. 



7.5.4 Density of Fisher zeros for the Ising model 

One trademark of F. Y. Wu's research is that he often looks at old problems 
and finds new life or new solution that others have not previously seen. A good 
example of a new look at an old problem is the density of Fisher zeros for the 
Ising model, which is the q = 2 Potts model. 

In 1964 Fisher pointed out that in the thermodynamic limit the zeros of 
the Ising partition function for a square lattice lie on two circles, now known 
as the Fisher circles, in the complex tanhii" plane, where K is the nearest- 
neighbor interaction. However, Fisher had not made the argument rigorous 
and, furthermore, no one has bothered to look into the distribution of the zeros 
on the circles, except at small angles which dictates the Ising critical behavior. 

Both of these two deficiencies have been rectified by W. T. Lu and F. Y. Wu. 
First, by considering zeros of the Fisher zeros for the Ising model on finite self- 
dual lattices, Lu and Wu ]190[ established rigorously that, indeed, the Fisher 

24 Here, we use the term circuit in the topological sense that solid triangles can be regarded 
as stars having three branches, each of which can be connected to other triangles and bonds. 
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zeros approach two circles in the thermodynamic limit. In a subsequent paper 
published in 2000 |202|| , they deduced the close-form expression for the density 
of the Fisher zeros for many regular two-dimensional lattices, thus completing 
the story of the Fisher zeros some 25 years after it was proposed! 



7.6 Duality relation for Potts correlation functions 
7.6.1 Correlation dualities 



Duality considerations are not often applied to correlation functions [133]. But 
F. Y. Wu initiated a new method for generating duality relations for correlation 
functions of the Potts model on planar graphs. In a pioneering paper [ 183 1 
he obtained duality relations for 2- and 3-point correlation functions, for spins 
residing on the boundary of a lattice. The consideration was soon extended to 



n-point correlations [187, 189 and to the case where the spins reside on 2 or 
more faces in the interior of the lattice [206]. A graph-theoretical formulation of 
the results in terms of rooted Tutte polynomials (se e sec tion ( 7.6.2[ ) below) was 
also given [186, 187 1. In addition, C. King and Wu |20f ] showed that generally 
it is linear combinations of correlation functions, not the individual correlations, 
that are related by dualities. 



7.6.2 Correlation function as rooted Tutte Polynomials 

As previously mentioned, the Potts partition function is also the Whitney- Tutte 
polynomial, or the Tutte polynomial for brief, considered in graph theory. In one 
further step, F. Y. Wu, C. King and W. L. Lu formulated the Potts correlation 
function as a rooted Tutte polynomial |l93fl . 

In graph theory a vertex is rooted if it is colored with a prescribed (fixed) 
color and a graph is rooted if it contains a rooted vertex. Interpreting the color 
of a site (vertex) as spin states, then the Potts correlation functions for which 
the spin states of given sites are fixed can naturally be formulated as rooted- 
Tutte polynomials. This is the basis of their graph-theoretical formulation of the 
Potts correlation function, from which dua lity relations of the Potts correlations 
become transparent and can be analyzed |l93||. 



7.7 Potts model and graph theory 

F. Y. Wu has written several review papers dedic ated to th e analysis of the Potts 
model from a graph theoretical viewpoint [^7], [71, 117, 122| . The most important 



one is the classic paper |37j written in collaboration with Baxter and Kelland 
alluded to above, which gives the graphical construction of the equivalence of 
the partition function of the Potts model with an ice-type model p7| . This 
derivation simplifies the algebraic method of Temperley and Lieb which is based 
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on the Temperley-Lieb algebra^]: 



if \i-j\>3, 



and applies to an arbitrary planar graph. This then opens the door for analyzing 
the triangular and honeycomb Potts models. 

Another important graphical analysis of a Potts model is the joint work with 
J. H. H. Perk [103, 104 on the non-intersecting string (NIS) model of Stroganov 
and Schultz, the close-packed loop model. The NIS model formulated by Perk 



and Wu in [ 103 1 turns out to be nothing but the bracket polynomial introduced 
by L. H. Kauffman in his state-model formulation of knot invariants. This fact 
offers a most natural approach to knot invariants from a statistical mechanical 



viewpoint [154] 



8 Other miscellaneous topics 

F. Y. Wu has worked on a diverse array of topics in mathematics and mathe- 
matical physics. In this section we present a random choice of topics that are 
not included above. 



8.1 Topics in graph theory 

F. Y. Wu is fond of graphs and has made many contributions to graph theory. 
A fine example is the aforementioned introduction of the rooted Tutte polyno- 
mial. Even when Wu does not obtain new results, he tries to provide simpler 
derivations, or find new consequences of known results. A good example is his 
work on random graphs |92| . 

Random graphs is a topic well-known to graph theorists after the work by 
Erdos and Renyi who introduced the problem in 1960. In the simplest formu- 
lation each pair of points of a set of N can be connected (by a bond, say) with 
a probability a/N, where a is a constant. One then asks questions such as 
what is the mean cluster size and the probability P(a) that the set becomes 
fragmented, namely not connected, etc. 

Using a Potts model formulation, Wu (9^[ reproduced the Erdos- Renyi result 
in just a few steps, showing that a transition occurs at a c = 1 and computed 
the critical exponents as well as the mean cluster size at criticality. This work 
has drawn considerable attention from graph theorists. Wu has also applied the 
result to evaluate the reliability probability of a communication network [ p2| . 

As another example of Wu's work in graph theory, one can mention his paper 
on the Temperley-Nagle identity for graph embeddings p9| where he provides 

25 A matrix representation of the C/j ,±i is for instance the q n X q n matrices with entries 

i- i/2 w ] =i;^ l ^j<°' J ) and g i/2 ^,^+i)n;=i' 5 (^.^)- 
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a simple derivation of the Temperley-Nagle identity 

= (i+^E^-i) ; (rb) u 

G L 

where G denotes section graphs of the original graph of v vertices and I lines, 
and L is a line set of G containing I lines covering v vertices. Using a weak- 
graph expansion he also deduced a sum-rule relation connecting the lattice con- 
stants of weak and strong embeddings. 

8.2 The vicious neighbor problem 

Consider N points randomly distributed in a bounded d-dimensional space. At 
a given instance, each point destroys his nearest neighbor (vicious neighbors) 
with a probability p. What is the probability Pn{p) that a given point will 
survive in the N — ► oo limit? 

The d = 2, p = 1 version of this problem was first posed by the Brandeis 
mathematician R. Abilock in American Mathematical Monthly in 1967, which 
remained unsolved for almost two decades. In 1986 the Omni magazine posted a 
prize for its correct solution, and R. Tao and Wu claimed the prize by publishing 



the solution for general d and p in 1987 |111 |. 

The idea of their solution is very simple. In d dimensions a given point can 
be killed by at most a finite number of other points. In two dimensions, 
for example, the number is n-i = 5 (theoretically a point can also be killed by 
6 other points but the phase space for that to happen has a zero measure). 
Therefore, one computes the volumes of the phase space for a point to be killed 
by 1, 2, • • •, rid neighbors, and uses the inclusion-exclusion principle to write the 
probability in question as an alternate series in p, whose highest power is n^. 
However, the evaluation of the volumes of the phase space is tedious requiring 
special techniques. 

For d = 1 the result is quite simple and one has Poa{p) = 1 — p + p 2 /2. For 
d = 2 the result is: 

Poo(p) = 1-0+0.316 3335p 2 -0.032 9390 /+0.000 6575/-0.000 OOlOp 5 (72) 

where the coefficient of each term is evaluated from integrals which can in prin- 
ciple be computed to any numerical accuracy. The coefficient of the last term, 
for example, is obtained by combining two 8-fold integrations. 

For p = 1 , @ yields Poo(l) = 0.284 051..., a solution which claimed the 
Omni prize. Tao and Wu also carried out Monte Carlo simulations to obtain 
the solution for d = 3, 4, 5. 

8.3 Counting partitions: from Potts to three-dimensional 
enumeration and beyond 



F. Y. Wu et al. [172 considered a directed lattice animal problem on the d- 



dimensional hypercubic lattice, and established its equivalence first with the 
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infinite-range Potts model and, in a second step, with the enumeration of (d— 1)- 
dimensional restricted partitions of an integer. The directed compact lattice 
animal problem was solved exactly in two and three dimensions, using known 
results in number theory. They found that the number of lattice animals of n 
sites grows as: 

exp(c-n (d - 1)/d ). 

Furthermore, the infinite-state Potts model solution leads to a conjectured lim- 
iting form for the generating function of restricted partitions for d > 3, which 
is a long-standing unsolved problem in number theory. 

Let us denote by A n (L\, L2, ■ ■ ■ Ld) the number of rt-site animals that can 
grow on an Li x L 2 x ■ • • x lattice. F. Y. Wu et al. showed that A n is 
precisely the number of (d — l)-dimensional restricted partitions of the integer 
n into non-negative parts to units of a hypercube of size L\ x L2 x ■ ■ • x Ld-i, 
with the size of each part at most Ld- 

Define the generating function 

L 1 -L 2 -L 3 ---L ( i 

G(Li, L%, L 3 , ■ ■ ■ Ld; t) = 1 + ^ A n (Li, L 2 , • ■ ■ Ld) ■ t n , 

n = l 

which is of interest in number theory. F. Y. Wu et al. showed that G is 
precisely the partition function of a Potts model on the d-dimensional lattice 
in the infinite-state limit, provided one identifies t with x d where x — (e — 
1)/ 'q x l d . This then connects the Potts model with the theory of partitions in 
number theory. 

For d — 2 the generating function corresponding to the square lattice 
(Li, L 2 ) readsQ: 

G(L u L 2 ;t) = (73) 

where (t) p — Yiq = i (1 — t q )- G(L\, L2) is a polynomial in t also known as the 
Gaussian polynomial or the "^-coefficient" . For d = 3 the generating function 
reads: 

G\L\, L-i, L 3 , t) - — — =-= , (74j 

[t\L 1 +L 2 ■ mL 2 + L 3 ■ [t\L 3 +Li 

where [i]x, denotes: 

L-l p 

[Al = II (*)p» L > 1, (t) p = J] (l-t«). (75) 

p=l q=l 



° Let us recall the classic analysis due to Rademacher which yields the celebrated Hardy- 
Ramanujan asymptotic result: A n — 1/(4 n ■ cxp(7r-^/ 2 n/3). One has the following 
asymptotic behavior for A n (L±, L2) when L\ = L2 : A n (L, L) ~ / (2 it n)) ■ 2 2 ^™. 
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For d — 4 G(Li, L 2 , L 3 , L4; t) is the generating function of restricted solid 
partitions of a positive integer into parts on a L\ x L 2 x i 3 cubic lattice, with 
each part no greater than L4. The evaluation of a closed- form expression for G 
in this case has remained an unsolved problem for almost a century. 

The expression which straightforwardly generalizes ((73|) and (|74|) would be: 

GstraightCii, L 2 , L 3 , L^-t) = t^tttt, with: (76) 

D G\t) 

N G(t) = {*}-£,!+ i 2 + i 3 + L 4 • {*}ii + i 2 ' {^}ii + i 3 " + Li 

-Dg(*) = {*}l 1 +L 2 + L 3 • {t}L 2 + L 3 + Li ■ {t}L 1 +L 3 + L 4 • +£ 3 + £4 

•{*}£! • • Wis ' {*}i< 



where: 



L-l 

= n Mp . x > 2 ( ?? ) 



But the explicit expression of G(2, 2, 2, L4; t) obtained by Major P. A. MacMa- 
hon in 1916 is 



G(2, 2, 2, L 4 ; t) - G straight (2, 2, 2, L 4 ; t) + G(2, 2, 2, L 4 ; t) (78) 



where 



G(2, 2, 2, L 4 ; i) = V 5l • .f L ;+ 8 ^ , with 



i=0 

.go =1, 91 = 2t 2 ■ (1 + t + t 2 + t 3 + < 4 ) + i 4 , 

i92 = i 5 ■ (1 +3i + At 2 + 8t 3 + 4t 4 +3i 5 + i 6 ) 
.93 - 2t 10 ■ (1 + t + t 2 + t 3 + t A ) + t 12 , g A = t 16 , 

and 



V < 2 + < + l / V(t) 8 (t) P _2- 

H. Y. Huang and F. Y. Wu [|l79| decided to look into the zeros of the 
generating function G(2, 2, 2, L 4 ; t) for various increasing values of i 4 . They 
found that the zeros are not exactly on the unit circle but seem to converge to 
the unit circle as L 4 increases. This indicates that a multiplicative correction 
G mu it(2, 2, 2, L 4 ; t) — G(Li, L 2 , L 3 , £4; t)/G str aight(-ki, ^2, ^3) ^4i*) would 
not have any simple Eulerian product form as in ( |7q) and j7^): 

G mult (2, 2, 2, L 4 ; t) = J] ( J fl ) Q " ( 79 ) 

n=l 

where a n are positive integers, since these product forms (|79| ) would necessarily 
yield zeros on the unit circle. H. Y. Huang and F. Y. Wu conjectured however, 
on the basis on their numerical results, that the zeros tends to be on the unit 
circle in the limit when any one of L%, L 2 , £3, L4 — * 00. 
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8.3.1 Directed percolation and random walk problems 

F. Y. Wu and H. E. Stanley J9(| have considered a directed percolation problem 
on square and triangular lattices in which the occupation probability is unity 
along one spatial direction. They formulated the problem as a random walk, and 
evaluated in closed-form the percolation probability, or the arriving probability 
of a walker. To this date this solution stands as the only exactly solved model 
of directed percolation. 



In another random walk problem Wu and H. Kunz 192] considered restricted 
random walks on graphs, which keep track of the number of immediate reversal 
steps, by using a transfer matrix formulation. A closed-form expression was 
obtained for the number of n-step walks with r immediate reversals for any 
graph. In the case of graphs of a uniform valence, they established a probabilistic 
meaning of the formulation, and deduced explicit expressions for the generating 
function in terms of the eigenvalues of the adjacency matrix. 



9 Knot theory 

The connection between knot theory and statistical mechanics was probably 
first discovered by Jones. His derivation of the V. Jones polynomial reflects 
the resemblance of the von Neumann algebra he uses with the Lieb-Temperley 
algebra occurring in the Potts model (see section (^^)). The direct connection 
came to light when L. Kaufman produced a simple derivation of the Jones 
polynomial using the very diagrammatic formulation of the non-intersecting 



string (NIS) model of J. H. H. Perk and F. Y. Wu |0§, [lOj. Soon thereafter 
Jones worked out a derivation of the Homfly polynomial using a vertex-model 
approach. The connection between knot theory and lattice statistical mechanics 
was further extended by Jones to include spin and IRF models. 

F. Y. Wu has written several papers on the connection between knot the- 



ory a nd s tatistical mechanics |150| , 151 , 154] including the comprehensive re- 



15C|. In hindsight, knot invariants arose naturally in statistical mechanics 



even before the connection with solvable models was discovered. In their joint 



paper 103 , for example, J. H. H. Perk and F. Y. Wu described a version of an 
NIS model which is precisely the bracket polynomial of L. Kauffman. Similarly, 
the g-color NIS model studied by J.H.H. Perk and C. Schultz is a q (q — 1) 
vertex model which generates the Homfly polynomial. Here we briefly describe 
the latter connection. 

The g-color NIS model has vertex weights (6 a bcd =1 if a = b = c = d and 
zero otherwise): 

w(a, b, c, d) = (80) 
(W(u) - S(u) - T(u)) ■ Sabcd + S(u) ■ S ab S cd + T(u) ■ S ac S bd , 

where: W{u) = sinh(u) = sinh(^ +it), S(u) = sinh(u), 

T{u) = sinh(?7 — u), q = e v + e~ v , 
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and the Homfly polynomial is a two variable knot invariant polynomial, discov- 
ered after Jones' work, by Freyd et al. The Homfly polynomial knot invariant 
has since been re-derived and analyzed by Jones using the Hecke algebra of 
the braid group. It can also be constructed from the Perk-Schultz NIS model. 
Actually the partition function Z(q, e v ) of the NIS model is a knot invariant 
related to the Homfly polynomial P(t, z): 

smh(<7?7) 

In the infinite rapidity limit this model leads to the Jones polynomial. The 
Boltzmann weight ( ^0| ) of the non-intersecting string model becomes: 

w{a, 6, c, d) = - e ±2r >6 a<b 6 c , d + e ±r >S b>d with: q = e" + e"" (81) 

The Jones polynomial V(t) is then obtained from the Homfly polynomial P(t, z) 
by taking z = \ft — 



FY. Wu discussed many knot invariants in his review [154]: the Alexander- 
Conway polynomial, the Jones polynomial, the Homfly polynomial, the Kauf- 
man polynomial and the Akutsu-Wadati polynomial, etc. The Alexander-Conway 
polynomial can be obtained from the Homfly polynomial by setting t = 1 in 
the Homfly polynomial P(t, z). The Akutsu-Wadati polynomial is an exam- 
ple of a new knot invariant derived from exactly solvable models in statistical 
mechanics. 

As our final example of F. Y. Wu's versatility, he and P. Pant and C. 



King [162 1, have obtained a new knot invariant using the exactly solvable chiral 
Potts model and a generalized Gaussian summation identity. Starting from a 
general formulation of link invariants using edge-interaction spin models, they 
establish the uniqueness of the invariant for self-dual models. They applied the 
formulation to the self-dual chiral Potts model, and obtain a link invariant in the 
form of a lattice sum defined by a matrix associated with the link diagram. A 
generalized Gaussian summation identity was then used to carry out this lattice 
sum, enabling them to cast the invariant into a tractable form. The resulting 
expression for the link invariant was characterized by roots of unity and does 
not appear to belong to the usual quantum group family of invariants. 



Finally, Pant and Wu [185| have derived a link invariant associated with the 
Izergin-Korepin model. 



10 Conclusion 

It would not be fair to summarize F. Y. Wu's contributions by a quick con- 
clusion such as: he wrote several important monographs on vertex models, on 
the Potts model and on knot theory, obtained many important results, in par- 
ticular the Lieb-Wu solution of the Hubbard model, the Fan-Wu free-fermion 
vertex model, the solution of the Baxter- Wu model, and many other results on 
dimers or free-fermion models, 3D dimers, d-dimensional free-fermion models, 
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Potts models, Ising and vertex models, using a large set of tools including ana- 
lytic calculations, expansions, series analysis, Monte-Carlo, with a particular 
emphasis on graph-theoretical methods. 

Most of the work of F. Y. Wu could be said to correspond to exact results 
in lattice statistical mechanics, or in mathematics, with particular emphasis on 
graph theory and enumerative combinatorics. We have tried to give here some 
hints as to F. Y. Wu's very large "graph" of concepts, results, tools, models, with 
many "intellectual loops" . We have not tried to provide an exhaustive descrip- 
tion of F. Y. Wu's contributions but, rather, only to provide a few comments 
on some of his results, emphasizing the fruitful cross-fertilizations between the 
various domains of mathematical physics and mathematics, and also to show 
the motivation and relevance of these results, tools, concepts and methods. 

Beyond a post-modern accountant's evaluation and from a research view- 
point, one must say that the important and numerous results F. Y. Wu has 
obtained arc not due to publish-or-perish productivity pressure, but, on the 
contrary, are the natural consequence of the pleasure of a scientist who loves to 
play with concepts and mathematical objects (dimers, graphs with particular 
boundary conditions, dualities, Potts models, series expansions with transmis- 
sivities, ...) and who has a strong desire to reach ambitious goals such as obtain- 
ing new results in three dimensions, new results for non-critical Potts models, 
or even for the Ising model in a magnetic field. 

A scientist does not become as productive as F. Y. Wu in response to external 
pressure but, on the contrary, only by forces being in harmony with his deep 
personal motivations. This is the only way to be as efficient and productive 
as F. Y. Wu and, as the famous French mathematician Jean Dieudonne once 
wrote, to work efficiently, "pour l'honneur de l'esprit humain" . 
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